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ABSTRACT 

Massive metal-poor stars might form massive stellar black holes (BHs), with mass 
25 < mBH/M0 < 80, via direct collapse. We derive the number of massive BHs (Nbh) 
that are expected to form per galaxy through this mechanism. Such massive BHs 
might power most of the observed ultra- luminous X-ray sources (ULXs) . We select a 
sample of 64 galaxies with X-ray coverage, measurements of the star formation rate 
(SFR) and of the metallicity. We find that Nbh correlates with the number of observed 
ULXs per galaxy (Nulx) in this sample. We discuss the dependence of our model on 
the SFR and on the metallicity. The SFR is found to be crucial, consistently with 
previous studies. The metallicity plays a role in our model, since a lower metallicity 
enhances the formation of massive BHs. Consistently with our model, the data indicate 
that there might be an anticorrelation between Nulx, normalized to the SFR, and 
the metallicity. A larger and more homogeneous sample of metallicity measurements 
is required, in order to confirm our results. 
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1 INTRODUCTION 

Ultra-luminous X-ray sources (ULXs, see Mushotzky 2004 
for a review, and references therein) are defined as non- 
nuclear point-like sources with isotropic X-ray luminosity 
ix Si 10^^ erg s~^. The mechanism that powers the ULXs 
is still unknown, although various scenarios have been pro- 
posed. ULXs could be associated with high-mass X-ray bi- 
naries (HMXBs) powered by stellar-mass black holes (BHs) 
with anisotropic X-ray emission (e.g. King et al. 2001) or 
with super-Eddington accretion rate/luminosity (e.g. Begel- 
man 2002; King & Pounds 2003; Socrates & Davis 2006; 
Poutanen et al. 2007) or with a combination of the two 
mechanisms (e.g. King 2008). ULXs could also be associated 
with HMXBs powered by intermediate- mass BHs (IMBHs), 
i.e. BHs with mass 100 M© < ttibh < 10^ M0 (see van der 
Marel 2004 for a review). IMBHs with mass larger than 
100 M0 may be required to explain the brightest ULXs (i.e. 
the < 4 ULXs with Lx ~ 10*^ erg s"^), those ULXs show- 
ing quasi-periodic oscillations (M82 X-1, see Strohmayer & 
Mushotzky 2003, and NGC 5408 X-1, see Strohmayer et al. 
2007) and some of those that are surrounded by isotropically 
ionized nebulae (PakuU & Mirioni 2002; Kaaret, Ward & 
Zezas 2004). However, IMBHs are not needed to explain the 
observational properties of most of the ULXs (e.g. Gongalves 



& Soria 2006; Stobbart, Roberts & Wilms 2006; Copper- 
wheat et al. 2007; Roberts 2007; Zampieri & Roberts 2009). 
Thus, the objects classified as ULXs might actually be an 
inhomogeneous sample, including sources of different nature. 

Most of ULXs are located in galaxies with a high star 
formation rate (SFR, e.g. Irwin, Bregman & Athey 2004), 
although a small fraction (10-20 per cent, especially in the 
low-luminosity tail) might be associated with an old popu- 
lation (e.g. Colbert et al. 2004; Brassington, Read & Pen- 
man 2005). The ULXs match the correlation between X-ray 
luminosity and SFR reported by various studies (Grimm, 
Gilfanov & Sunyaev 2003; Ranalli, Comastri & Setti 2003; 
Gilfanov, Grimm & Sunyaev 2004a,b,c; Kaaret & Alonso- 
Herrero 2008). Furthermore, the same studies indicate that 
the luminosity function of ULXs is the direct extension of 
that of HMXBs. Recent papers suggest a correlation be- 
tween ULXs and low-metallicity environments, and propose 
that this may be connected with the influence of metallicity 
on the evolution of massive stars (PakuU & Mirioni 2002; 
Zampieri et al. 2004; Soria et al. 2005; Swartz, Soria & Ten- 
nant 2008). This scenario has been explored in detail by 
Mapelli, Colpi & Zampieri (2009, hereafter Paper I) and by 
Zampieri & Roberts (2009), highlighting that a large frac- 
tion of ULXs may actually host massive (~ 30 — 80 Mq) 
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stellar BHs formed in a low-metallicity environment. In fact, 
low-metallicity massive ( ;J, 40 Mq) stars lose only a small 
fraction of their mass due to stellar winds (Maeder 1992, 
hereafter M92; Heger & Woosley 2002, hereafter HW02; 
Heger et al. 2003, hereafter H03; Belczynski et al. 2010, 
hereafter BIO) and can directly collapse (Fryer 1999; BIO) 
into massive BHs (25 M© < ttibh < 80 M©). Such massive 
BHs can power most of the known ULXs without requiring 
super-Eddington accretion or anisotropic emission. Further- 
more, their formation mechanism can explain the correlation 
between ULXs and SFR, and the fact that ULXs are pref- 
erentially found in low-metallicity regions. In this Paper, we 
extend to a larger sample of galaxies the analysis reported 
in Paper I, and we study the formation of massive BHs from 
the direct collapse of massive metal-poor stars and their 
possible connection with ULXs. In particular, we show that 
there is a correlation between the number of massive BHs 
formed in this scenario and the number of observed ULXs 
per galaxy. 



2 METHOD 

The aim of this Paper is to compute the number of massive 
BHs that are expected to form in a galaxy (Nbh) through 
the direct collapse of massive metal-poor stars. We will then 
compare such number with the observed number of ULXs 
(NuLx) in the same galaxy and with other observational 
quantities (e.g. the SFR, the metallicity, etc.). In this Sec- 
tion, we start describing the procedure adopted for deriving 
Nbh and then present the properties of the galaxy sample 
used for the comparison. 



2.1 The number of massive BHs 

According to numerical models (Fryer 1999; HW02; H03), a 
star that, at the end of its life, has a final mass mfln > 40 M0 
is expected to directly collapse into a BH. In this case, the 
mass of the remnant BH is likely more than half of the final 
mass of the progenitor star, as relatively small mass ejection 
is expected in the direct collapse. Thus, stars that at the 
end of their lives have mfln > 40 Mq are likely to produce 
massive BHs (BIO). 

The final masses of the stars strongly depend on their 
metallicity. Massive stars with metallicity close to solar can- 
not have final masses larger than mfln ~ 10 — 15 Mq, even 
if their initial mass was very large, as they are expected to 
lose a lot of mass due to stellar winds (H03). Instead, mas- 
sive stars with lower metallicity are less affected by stellar 
winds, and retain a larger fraction of their initial mass. If its 
metallicity is sufficiently low, a star can have a final mass 
Tifln > 40 Mq and can directly collapse into a massive BH 
with a mass 25 M© < thbh < 80 M© (HW02; H03; BIO). 

Thus, in order to compute the number of massive BHs 
Nbh, we first need to know for which metallicities massive 
stars can directly collapse into BHs. Evolutionary tracks 
of massive stars have been computed under different as- 
sumptions concerning the mass-loss rate (M92; Bressan et al 
1993; Fagotto et al 1994a,b; Eldridge & Tout 2004; Eldridge 
& Vink 2006; BIO), and thus predict different final stellar 
masses (differences are up to a factor of 2). In Paper I, we 
used the results from M92, that can be regarded as upper 



limits. In this Paper, we base our calculations on the results 
of two previous studies: one relies on the PAD OVA tracks 
(as described in Portinari, Chiosi & Bressan 1998, hereafter 
referred to as P98) and the other is presented by BIO. P98 
assume that all massive stars explode as supernovae, and 
do not consider the direct collapse scenario. Here, we follow 
the more recent studies by Fryer (1999) and by H03, indicat- 
ing that, above 40 Mq, the stars collapse to a massive BH. 
Thus, in the case of P98, we take the final stellar masses for 
various metallicities listed in their table 1, and we assume 
that all stars with men > 40 M© end up into a massive BH^ . 

Instead, BIO already include in their code the hypoth- 
esis that massive stars directly collapse into massive BHs. 
Thus, when we apply the results from BIO to our model, we 
consider the stars that end into BHs with rriBH ^ 25 Mq 
as progenitors of massive BHs (see fig. 1 of BIO). This 
condition is approximately equivalent^ to requiring that 
mfln >40M©. 

Given these assumptions, we can derive the expected 
number of massive BHs per galaxy (Nbh) as a function of 
the star formation rate (SFR) and of the metallicity Z (see 
Paper I): 



Nbh(SFR, Z) = A(SFR) 



' dm. 



(1) 



!{Z) 



where m,max is the maximum stellar mass (we assume 
m,max = 120 M©) and a is the slope of the initial mass func- 
tion (IMF), mprog(^) is the minimum initial stellar mass 
(i.e. the mass at zero-age main sequence) for which a star 
is the progenitor of a massive BH. As we discussed above, 
rripiogiZ) strongly depends on the metallicity. In our calcu- 
lations, we assume m,prog(^) to be the initial stellar mass for 
which mfln > 40 M© and msn ~ 25 M© , when adopting the 
models by P98 and by BIO, respectively. Finally, A(SFR), 
the normalization constant in equation (1), can be estimated 
as 



A(SFR) 



m^ ^ dm 



(2) 



where rrimin is the minimum stellar mass (we assume m-min ~ 
0.08 M©), SFR(t) is the star formation rate as a function of 
time and tspR is the duration of the star formation. 



^ An important difference from the original assumptions made by 
P98 concerns the fate of very massive stars, with Hydrogen ex- 
hausted core at the end of the hydrostatic evolution, Mjjoi larger 
than 40 M©. For Mhc >40 M©, P98 adopted the pair instabil- 
ity (Woosley & Weaver 1986) model with the following fates: for 
Mho in the range 40 to 60 M© the star undergoes a pulsational 
instability leading to a final iron core collapse; in the range Mhb 
~60 to ^100 M© the star was supposed to undergo a complete 
disruption and finally, for even larger cores, a direct collapse to a 
BH. 

■^ We assume rriBH ^ 25 M© for BIO, because the models (e.g 
HW02; H03) show that in the direct collapse slightly more than 
half of the final mass of the star goes into the BH. As the min- 
imum final stellar mass for direct collapse is rragn ^ 40 M© , 
"iBH ;CL 20 — 30 M© is a reasonable choice. On the basis of this 
argument, threshold masses anywhere in between 20 and 30 M0 
would be fine. We choose rriBH ~ 25 M© as a reasonable fidu- 
cial value. Furthermore, BHs with msH X. 25 M© are above the 
observed mass-range for stellar BHs (3 — 20M©, Orosz 2003). 
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However, we are not interested in all the massive BHs, 
but only in those that could acquire massive stellar com- 
panions and power observable ULXs. Thus, in the case of 
interest, equation (2) might be written as 

SFR tco 



A(SFR) 



J rri 



m^~°' dm 



(3) 



where SFR is the current star formation rate and tco is the 
characteristic lifetime of the companion of the massive BH. 
In this Paper, we adopt a constant value ico = 10^ yr. Such 
value is the lifetime of a ~ 15 M© star. In fact, according to 
theoretical models, massive BHs can power persistent ULXs 
only when their companion is sufficiently massive (> 15 M©, 
Patruno et al. 2005). If the companion star is relatively low- 
mass (2 < m/M0 < 10), then the system can power only a 
transient ULX, with an outburst phase of a few days every 
few months'^ (Portegies Zwart, Dewi & Maccarone 2004). 
The probability of observing such transient ULXs is a factor 
of ^ lO'^ lower than the probability of observing persistent 
ULXs. Thus, the value of Nbh derived assuming ico ~ 10^ yr 
is not the total number of massive BHs present in a galaxy, 
but the number of massive BHs that could easily find a 
massive companion to accrete from. 

Fig. 1 shows Nbh, normalized to the SFR, as a function 
of Z, for P98 and for BIO. In this Fig., we adopt two differ- 
ent IMFs: (i) the Salpeter IMF (Salpeter 1955), for which 
Q = 2.35, and (ii) the Kroupa IMF, for which a = 1.3 if 
m < 0.5 M0 and a = 2.3 for larger masses (Kroupa 2001). 
Fig. 1 indicates that Nbh has approximately the same trend 
when adopting the Salpeter or the Kroupa IMF, although 
the normalization is different by a factor of ~ 2. In the fol- 
lowing, we will consider, as reference, only the Kroupa IMF. 
As shown in Appendix A, where a more detailed comparison 
is reported, the results can be easily scaled to the Salpeter 
IMF. Fig. 1 shows also that the trend obtained adopting 
P98 is quite similar to the one obtained adopting BIO. The 
main differences are in the normalization, as Nbh is a factor 
of ~ 1.3 — 2.0 lower for P98 than for BIO, and in the cut-off, 
as Nbh goes to zero for Z ~ 0.4 Zq and 0.5 Zq for P98 and 
BIO, respectively. 

Once Nbh is known, we can estimate the upper limit 
(ebh) of the fraction of massive BHs that power ULXs in 
a given galaxy at present, assuming that all the observed 
ULXs in this galaxy, corrected for the contamination (Nulx, 
see next Section), are powered by a massive BH: 

Nulx 



ebh ~ 



Ne 



(4) 



2.2 The sample of galaxies 

In order to compare the value of Nbh computed from equa- 
tion (1) with the observed number of ULXs per galaxy 
(Nulx), we selected a sample of galaxies (Table 1 and refer- 
ences in Appendix B) that satisfy the following criteria, (i) 



■^ However, the Galactic transient source GRS 1915 + 105, which 
is associated to a low-mass X-ray binary system including a quite 
massive BH (mBH = 14.0 ± 4.4 M©, Harlaftis & Greiner 2004), 
has a much longer outburst phase (> 17 years, see e.g. Castro- 
Tirado, Brandt & Lund 1992; Fender & BeUoni 2004; Deegan, 
Combet & Wynn 2009), in partial contradiction with this model. 
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Figure 1. Number of expected massive BHs per galaxy, nor- 
malized to the SFR, as a function of Z. Sohd line (red on the 
web): BIO with the Kroupa IMF; dashed line (green on the web): 
BIO with the Salpeter IMF; dot-dashed Une (black): P98 with the 
Kroupa IMF; dotted line (blue on the web): P98 with the Salpeter 
IMF. 



There exist X-ray observation(s) deep enough to detect the 
presence of ULXs. (ii) The selected galaxies have at least 
one measurement of the SFR. (iii) They have a sufficiently 
accurate measurement of the metallicity Z (in most cases 
even of the metallicity gradient), (iv) They are not elliptical 
galaxies. 

The last criterion indicates that the galaxies in Table 1 
are spiral, irregular or peculiar galaxies, but not elliptical 
galaxies. We decided not to include any elliptical galaxy in 
our sample, as ULXs in elliptical galaxies show properties 
that are quite different from those of ULXs in other galax- 
ies, which may indicate a different origin for them (Swartz et 
al. 2004; Liu, Bregman & Irwin 2006; Winter, Mushotzky & 
Reynolds 2006). Furthermore, we remind that in early- type 
galaxies there is neither evidence of the presence of core- 
collapse supernovae nor of a significant fraction of massive 
star population, apart from small rejuvenation episodes (e.g. 
NGC 205, Bertola et al. 1995). ULXs in elliptical galaxies 
might be connected to the minor (but non- negligible) frac- 
tion of ULXs that is claimed (e.g. Colbert et al. 2004; Brass- 
ington et al. 2005) to be associated with old stellar popula- 
tions (hereafter, we refer to them as Population II ULXs). 
Since old populations are present also in spiral galaxies, our 
sample may be contaminated by Population II ULXs. It is 
quite difficult to estimate their number, but it is reasonable 
to neglect their contribution in our sample (see Appendix B 
for details). 

No complete catalogue of galaxies with and without 
ULXs exists. The most known available catalogues are Liu & 
Mirabel (2005; hereafter LM05) and Liu & Bregman (2005; 
hereafter LB05). LM05 is a catalogue of ULXs reported in 
the literature until April 2004. It does not include galaxies 
without ULXs, and it is quite non-uniform, as it reports ob- 
servations with various instruments. LB05 is a catalogue of 
all the 313 nearby galaxies that have been observed with the 
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Table 1. Properties of the galaxies in our sample. 



Galaxy 


SFR (Mq yr-i; 


) Z (Zq) 


NuLX.raw ^ 


D (Mpc) b 


Q" 


NULX'' 


The Cartwheel 


20 


0.14 


19 


124 





19.00+^-g 


NGC 253 


4.0 


0.24 


3 


3.0 


n 02+°'°'* 

"■"^-0.01 


2.9812-91 


NGC 300 


0.14 


0.19 





2.02 





o-oolJ:gg 


NGC 598 (M 33) 


1.1 


0.32 


1 


0.92 


n on+°°° 


1.00+2-29 


NGC 628 (M 74) 


2.2 


0.27 


2 


8.5 


-R-S§ 


1.85l2-'63 


NGC 1058 


0.27 


0.26 


1 


9.85 


n nq+°°s 


0-9llo;83 


NGC 1073 


1.2 


0.34 


2 


15.2 





2.0012-63 


NGC 1291 


0.94 


0.06 


3 


8.6 


16+°-" 
"•-■-"-0.07 


2 84+2-91 


NGC 1313 


1.4 


0.1 


2 


3.7 





2.00l2-g 


NGC 1365 


7.2 


0.20 


13 


17.4 


0-821°:?° 


12-18ltl^ 


IC342 


0.48 


0.19 


2 


3.9 


o.o2i«:°t 


l-SSli'i 


NGC 1566 


3.2 


0.33 


4 


13.4 


i-5ot«:l« 


2.50l?;i6 


NGC 1705 


0.087 


0.19 





6.0 





o^oolJ:^^ 


NGC 2366 


0.075 


0.10 





2.7 





o^ooljf,?, 


NGC 2403 


0.4 


0.22 


1 


3.1 





i^oolg;f3 


NGC 2442 


4.6 


0.45 


2 


17.1 


-1 OO + O.IO 

^••^•^-0.31 


0.67l2-f, 


Holmberg II (Arp 268) 


0.1 


0.1 


1 


4.5 





l-OOlg.fs 


NGC 2903 


1.9 


0.28 


2 


7.4 


n nq+°°'^ 


1 qi+2.63 
1-91_1.29 


NGC 3031 (M 81) 


3.3 


0.29 


2 


3.4 


n n3+°°5 


1.9712-63 


NGC 3049 


0.57 


0.60 





19.9 





o-oolJ:gg 


PGC 30819 (IC 2574, UGC 05666) 


0.057 


0.17 





3.5 





o^ooljf,^ 


NGC 3310 (Arp 217) 


2.2 


0.22 


3 


18.7 


1fi+°-03 
"•^"-0.03 


2.841-.S 


NGC 3395/3396 (Arp 270) 


4.7 


0.21 


7 


24.6 





7.00l^-i« 


PGC 35286 (UGC 06456) 


0.017 


0.062 





1.4 





o^oolJ:^^ 


PGC 35684 (UGC 06541, Mkn 178) 


0.01 


0.091 





4.0 





o^oolJ:f,f, 


NGC 3738 (Arp 234) 


0.038 


0.20 





4.3 





o^oolJ:gg 


NGC 3972 


0.30 


0.35 





17.0 





o^oolJ:^^ 


Antennae (NGC 4038/4039, Arp 244) 


7.1 


0.04 


15 


21.5 





15.001*:^^ 


NGC 4144 


0.05 


0.20 





5.7 





o^oolJ:^^ 


NGC 4214 


0.13 


0.21 





3.5 





o^oolJ:f,f, 


NGC 4236 


0.12 


0.22 





2.2 





o^oolJ:f,f, 


NGC 4248 


0.018 


0.26 





7.3 





O^OO-o'.oo 


NGC 4254 (M 99) 


4.0 


0.26 


3 


16.8 


1 QO+O.ll 

j^-^y_o.i3 


1 61+2-91 


NGC 4258 (M 106) 


2.5 


0.32 


3 


7.7 


0.58+g-« 


2 42+2.92 


NGC 4303 (M 61) 


5.8 


0.20 


6 


15.2 


1 Q1+0.58 


4 09+3-'^° 
4.uy_2 45 


NGC 4321 (M 100) 


4.8 


0.28 


7 


14.1 


1 SS+0-65 
t-»»-0.49 


c 12+3-79 


NGC 4395 


0.12 


0.22 


1 


3.6 


n n4+°-°'^ 

"•"^-0.03 


n Qf.+2.29 


NGC 4449 


0.28 


0.25 





3.0 





o^oolJ:f,f, 


NGC 4485/4490 (Arp 269) 


4.5 


0.23 


5 


8.55 


16+°-" 
"•-■-"-0.06 


4.84l^-^« 


NGC 4501 (M 88)° 


3.5 


0.49 


2 


16.8 


"•^•^-0.09 


1 57+2-63 
— 1.29 


NGC 4559 


1.3 


0.19 


2 


5.8 





2.0012-63 


NGC 4631 (Arp 281) 


2.80 


0.22 


2 


6.9 


n i6+°-" 


^■84_i.3o 


NGC 4651 (Arp 189) 


1.40 


0.23 


1 


16.8 


17+0-04 
"-^'-0.03 


r, oq+2.29 


NGC 4656 


0.95 


0.075 


1 


7.2 


05+0-04 
"•"0-0.02 


0^95l2:2«3 


The Mice (NGC 4676, Arp 242) 


6.0 


0.3 


6 


90.4 





6.00+^-^^ 


NGC 4736 (M 94) 


1.1 


0.20 


1 


5.2 


15+°-^'^ 

"--"-0-0.08 


f. ot:+2:29 
^■8&-0.85 


NGC 4861 (Arp 266) 


0.59 


0.13 


2 


17.8 


11+0-03 
"-^^-0.02 


1.89l2-« 


PGC 45561 (UGC 08231) 


0.29 


0.40 





33.0 





o-ool,V.f,f, 


NGC 5033 


1.50 


0.069 


2 


15.2 


omtlil 


1 01+2.63 
1-^1-1.30 


NGC 5055 (M 63) 


1.70 


0.30 


2 


8.5 


o.i8lg:J? 


1.8212-63 


NGC 5194/5195 (M 51, Arp 85) 


13.0 


0.39 


9 


7.7 


11+0-08 
"•^^-0.05 


8.891^-10 


NGC 5236 (M 83) 


2.6 


0.36 


1 


4.7 


n 09+0-12 

^••Jy-o.os 


O-Sllo'.M 


NGC 5238 (Mkn 1479) 


0.013 


0.26 





4.9 





o-oolJ:f,f, 


NGC 5408 


0.09 


0.11 


1 


4.85 





i-oolg;g 


NGC 5457 (M 101, Arp 26) 


3.1 


0.17 


7 


6.9 


n 40+0-34 


6-6ol^:S 



Table 1 - continued. 
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Table 1 - continued. 



Galaxy 



SFR (M0 yr-l) Z (Zq) 



D (Mpc) b 



Circinus 


1.5 


0.1 


4 


NGC 6946 (Arp 29) 


3.56 


0.283 


3 


PGC 68618 (IC 5201) 


1.7 


0.260 


1 


NGC 7714/7715 (Arp 284, Mkii 538) 


7.2 


0.2 


9 


NGC 7742 


1.27 


0.245 


2 


Milky Way 


0.25 


0.306 





IC 10 


7.14 X 10-2 


0.22 





Large Magellanic Cloud (LMC) 


0.25 


0.27 





Small Magellanic Cloud (SMC) 


0.15 


0.129 






4.2 

5.5 

11.1 

36.6 

22.2 

0.70 

0.051 

0.064 



0.00+0 

O.OSto' 

0.49t°- 

0.84t° 

0.32to, 











4.00 
2.92 
0.51 
8.16 
1.68 
0.00 
0.00 
O.OO 
0.00 



+3.15 
1.91 

+2.91 
1.63 

+2.29 
0.51 

+4.11 
2.94 

+2.63 
1.29 
1.83 
0.00 
1.83 
0.00 
1.83 
0.00 
1.83 
0.00 



Notes. The references for the values of SFR, Z and Nulx raw reported in this Table are given in Appendix B. '^ Nulx raw is the 

observed number of ULXs per galaxy, before the subtraction of contaminating sources. '^ D is the distance in Mpc (from the Liu & 

Bregman 2005 catalogue; if the galaxy is not in such catalogue, from the NASA Extragalactic Database). '^ Number of expected 

contaminating objects, estimated as described in Sections 2.2.1 and B4 (Appendix B). When Q is equal to zero (without uncertainties), 

either the contamination was already subtracted in Nulx raw (i-e. in the paper with the original data) and Nulx raw = Nulx> or 

Nulx, raw = 0. '^ Nulx is the observed number of ULXs per galaxy, after the subtraction of the estimated number of contaminating 

sources. The errors come from the combination of the Poissonian uncertainties upon Nulx raw (according to the treatment in Gehrels 

1986), and on Q (with uncertainties from the Hasinger et al. 1998 log(N) — log(S) relation). ° NGC 4501 is the only galaxy with 

Nulx > and with Z > 0.47 2©. For this galaxy Nbh = in both P98 and BIO models (see Fig. 1 and Section 3.2), although 

Nbh + Icr > 0. 



RO S AT Yii^ Resolution Imager (HRI) with detection lim- 
its that would allow the detection of ULXs. Thus, the LB05 
catalogue is less biased than LM05, as it includes galaxies 
with and without ULXs, and it is more homogeneous, as 
all the data reported come from the same instrument. On 
the other hand, even LB05 suffers from some observation 
bias (i.e., galaxies were originally chosen to be observed for 
a myriad different reasons), and the instrumental properties 
(especially sensitivity and resolution) of ROSAT are quite 
worse than those of Chandra and XMM-Newton, whose ob- 
servations are included in LM05. Furthermore, some galax- 
ies have been observed after 2004 and are not included in 
these catalogues, but there is no reason not to take them 
into account. 

Given this situation, our sample is inevitably affected by 
various biases (e.g. galaxies without ULXs are likely under- 
represented) and it is not homogeneous. It is mainly based 
upon the LB05 catalogue, which includes 52 of the 64 galax- 
ies listed in Table 1. The remaining 12 galaxies can be di- 
vided into 5 Local Group galaxies (NGC 598, the Milky Way, 
IC 10, and the two Magellanic Clouds), and 7 non-Local 
Group galaxies (the Cartwheel, the Antennae, the Mice, 
NGC 628, NGC 1058, NGC 5408, and Circinus). The 5 Local 
Group galaxies have low SFR (the average is 0.36 M© yr^^, 
compared to 2.06 M© yr^^ for the 52 LB05 galaxies) and a 
total of only one ULX (i.e. an average of 0.2 ULXs/galaxy, 
compared to 2.29 ULXs/galaxy in the 52 LB05 galaxies; here 
we are neglecting contamination); whereas the 7 non-local 
group galaxies have high SFR (average of 5.31 M© yr"^) and 
many ULXs (average of 6.86 ULXs/galaxy). 

Each of the two non-LB05 sub-groups is clearly different 
from the sample of LB05 galaxies. However, we still include 
both of them in our sample, because they allow us to explore 
regimes that are not well represented by the 52 galaxies in 
LB05: the Local Group galaxies enlarge the number of galax- 
ies with no ULXs, that are likely under-represented because 
of observation and publication bias; the non-Local Group 



galaxies provide us with examples (such as the Cartwheel 
and the Antennae) of relatively rare objects with high SFR. 
We note that, even if a galaxy is in the LB05 catalogue, 
we often use X-ray measurements coming from more recent 
papers where dedicated analysis have been published: Ta- 
ble 1 includes 27 galaxies for which the raw number of ULXs 
(Nulx, raw, i.e. the observed number of ULXs per galaxy, be- 
fore the subtraction of contaminating sources) is taken from 
LB05, 24 galaxies for which it is taken from LM05, and 13 
galaxies for which Nulx, raw is taken from other papers. 



2.2.1 Short description of Table 1 

A full account of the origin of all the data reported in Table 1 
will be given in Appendix B. Here we only provide a basic 
description of the most important quantities. 

The SFRs come from either ultra-violet. Ha, far- 
infrared, radio measurements, or an average of these dif- 
ferent measurements, depending on the data available for 
each galaxy. Galaxies whose SFR is highly uncertain (a fac- 
tor of 10 or more), due, e.g., to the presence of an active 
galactic nucleus (AGN), are not considered in our sample 
(e.g. NGC 1068, NGC 3690). Concerning metallicity mea- 
surements, their existence is the most restrictive among our 
selection criteria, as metallicities are often unavailable. For 
most galaxies in the sample, we use measurements derived 
from line intensities in HII regions, translated into abun- 
dances with the empirical relation of Pilyugin (2001a; here- 
after POl) or Pilyugin & Thuan (2005; hereafter PT05). If 
more than one HII region was measured within a galaxy, it is 
generally possible to derive the metallicity gradient; in such 
cases, we use the metallicity at the mean distance of ULXs 
from the centre of their host galaxy (0.73 i?25, where _R25 is 
the radius of the 25th magnitude isophote; see Appendix B 
and Liu et al. 2006). When the spectra of HII regions are 
unavailable, we use X-ray measurements, although they are 
much less accurate. 

The number of expected contaminating sources (Q) 
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Table 2. Values of Nbh and esH for the galaxies listed in Table 1, assuming the models by P98 (first and second column) and by BIO 
(third and fourth column). We adopt a Kroupa IMF. 



Galaxy 



Nbh (P98) <:bh/10-* (P98) Nbh (BIO) eBn/lQ-^ (BIO) 



n+21700 
^-21600 

6070+"^^*'' 
UU/U_324Q 

245+125 
1280_;^280 

2860li«5™ 

077+198 
•^' '-203 
12Q0+^^'' 

ZiOU_j^3gQ 

3460+1760 
12600+®''°° 

356011^,^° 
1531^? 

"^' -339 

r,+5270 

247ti25 

2440ll|« 

43001^3^° 

Q + O 

3510llf,° 
74q0+'l°''° 

48+24 

25+" 

Z0_j^3 

63133 

01Q+174 
22200+1"°° 

ZZZUU_jj200 

"^-45 
910 + 115 

^^•-"-iio 
191+100 

25 +1* 
^0 -14 

54801^880 
290011610 

g7QQ + 5050 

^"-"J-5150 
g24Q+3420 

197tlgl 

4081^1^ 

6870l360g 

n+3180 

'-'-0 

9q7n+1230 
ZO(U_j^260 
4450+2400 

2080llJ3g 
9550+1290 
698013870 

1930+9*° 

1 300+'^70 

Q+377 

405Q+2050 

2040+1"° 
ZU*U_2040 

r, + 17300 

'-'-0 

r,+3580 

'-'-0 

-,0+10 

1^-10 

OOQ + 113 
^.^O 1 1 Q 

5890"1^??° 



-35200 
-35400 
11200+^700 

438li^ 

2510+1330 
Z01U_j3gQ 

5440+288° 
709^374 

2620+139° 

oggQ + 1860 
5100+2580 
20900+11°°° 

1500lJ«(! 

678013730 

274+139 

979+138 
^'^-136 

1160+*^°° 
ll°^-600 

5930tg70 

364ll«* 
4640+2450 

sisn+4330 

Q + O 

^-0 

-,07+95 

638013^*° 
1360017000 

681^4 
3811^ 

ll^lgg 

628133^* 

30000liJ°°° 

151176 

4041^?* 
3481176 

481^^ 

10700lfg°° 

569013010 

16800l«900 

11900+6300 

359+185 

750+400 

1 2600+6500 

n+6710 
^-0 

4240+2150 
'iZ'iU_22oo 

809011160 

3970+2010 

'^'790+1870 
■5I^'-'_1870 
13700+7500 
id^UU_72gg 

3310ll?J° 
2080ll°« 
5311^95 
5920+2970 

3990+211° 

0»»U_2ggg 

9qqr,ri+13600 
ZOOUU_233gg 

4000+2660 
'ia»u_3g2g 

35+18 

090+166 
32»_ig4 

gg^Q+4870 

»"J^'-'_4860 



The Cartwheel 
NGC 253 
NGC 300 
NGC 598 
NGC 628 
NGC 1058 
NGC 1073 
NGC 1291 
NGC 1313 
NGC 1365 
IC 342 
NGC 1566 
NGC 1705 
NGC 2366 
NGC 2403 
NGC 2442 
Holmberg II 
NGC 2903 
NGC 3031 
NGC 3049 
PGC 30819 
NGC 3310 
NGC 3395/3396 
PGC 35286 
PGG 35684 
NGC 3738 
NGC 3972 
Antennae 
NGC 4144 
NGC 4214 
NGC 4236 
NGC 4248 
NGC 4254 
NGC 4258 
NGC 4303 
NGC 4321 
NGC 4395 
NGC 4449 
NGC 4485/4490 
NGC 4501 
NGC 4559 
NGC 4631 
NGC 4651 
NGC 4656 
The Mice 
NGC 4736 
NGC 4861 
PGG 45561 
NGC 5033 
NGC 5055 
NGC 5194/5195 
NGC 5236 
NGC 5238 
NGC 5408 
NGC 5457 



419001 



4.5 



+2.7 
-2.6 

-5.4 



4 9^ 

0+74-7 

^-0 

7 0+18.4 

'•^-7.8 

6 7+9-9 

— 5 9 

24 1+62-1 
^^•1-24.1 
-15 5+22.2 
10-''-15.5 

10.6111^9 

5 8+8-1 

6 5+*-7 
"■0-4.2 

+33.7 

19.7 

7 n+10.3 

'■^-7.0 

n+119.6 



23.5"* 



-0 
+98.9 



+36.9 



15.7__^ 

> 1.3 

40 4+95.00 

D + 11.6 



7.8^ 



7.0 

6.6 







.r 



+ 166.4 


+9.4 
6.4 
9 3+7-1 
0+381.2 
^-0 
0+732.0 
^-0 
O+290.5 
^-0 
0+59. 03 
'-'-0 

6-7l^:^ 

n+217.9 
'-'-0 

0+85.9 
^-0 
+95.8 


+732.0 


+5.9 
2.9 
11.2 
8.3 
+4.5 
3.6 
o 9+7.8 
8-^-6.4 

48^7+119.2 







2.9 

8.3 

4.2 




7.1 



48.7 
+44.8 


+6.2 


>4.1 
8 4+11-9 

4 1+6-3 
*-l-3.7 
4 0+11-2 

7+9.2 
"^■'-3.7 

8.6ll^ 
4.4lf4i 
14 6+'2l-6 
-,+48.5 
-0 

1+6.9 

-3.4 

3+13.8 

-8.9 



0: 

3.4^ 
8.9J 

> 5.2 

> 2.5 

0-|-1017.0_g 

44 8+105.2 

,+8.6 

-7.4 



11.2"' 



70200: 



9 7+1.6 
^■'-1.5 
7+2.9 
-2.0 
-,+41.8 
-0 

-,+9.4 

-4.0 

-+5.2 

-3.0 

s+33.0 

-12.8 

™ + 10.8 

-6.5 
:i+8.8 
-6.0 
o q+5.5 
o.»_3 2 
a+3.8 
-3.4 
-18.7 
-11.0 
7+5.4 
-3.7 
-,+66.8 

-0 
-,+67.3 
-0 
.,-1-20.2 

i.4 
,+4.9 
-1.1 
-+64.4 
" -26.6 
,+6.1 
-3.6 



2.7^ 

Oj 

4.0! 

3.5l 

12. 8^ 

7.6l 

7.9J 

3.9l 

5.8l 

13.2! 

3.7:* 

0^ 

oj 

8.6! 

l.lj 

27.5^ 

4.1-' 



2.4 



+3.5 
-2.0 



,+97.9 

-0 

+5.1 
-3.5 

+3.8 

3.2 

,+269.1 


,+481.6 



■160.5 


,+29.1 







4.4 

5.1 











5.0 



1+3.0 
-2.8 
121.2 



0-0 
n+45.3 

^-0 

1+52.6 







-0 







,+381.2 



1 5+3.0 
l-°-1.5 
40+5.7 

+2.6 

2.0 

+4.1 

-3.3 

26.7+65-3 



2.4! 
4.3"' 



0! 
3.9! 
> 1 

4.7j 
2.3"' 



-27.5 
r,+24.1 
-0 

n+3.3 
-2.7 

> 1.9 

7+6.6 

-3.9 
0+3.5 

-2.0 
9 1+5.9 
^•1-2.1 
2 6+6-3 
4 43.6 

2 6+7-0 
^•"-2.6 

9 1+13-* 

Q+34.5 
^-0 
2 3+''-7 

4 6+7-0 

Q+2.8 
3-0-3.8 

1 0+4.7 
1-8-1.8 
Q+522.9 



-0 

30.5 
6.9 



71.5 
29.6 

,+5.2 

-4.4 



Table 2 - continued. 
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Table 2 - continued. 



Galaxy 



Nbh (P98) eBH/lO-* (P98) 



N 



BH 



(BIO) eBH/10-4 (BIO) 



Circinus 


07-10+1880 


10.8li"5i 


NGC 6946 


.f,of>+2540 

4DdU_2ggQ 


6-3lI:^ 


PGC 68618 


233011^20 


2 2+1'''-'- 


NGC 7714/7715 


12600+««» 


q 7+6-2 
^■'-5.9 


NGC 7742 


loqn+lOOO 

l»3U_ggQ 


9.2il«i3 


Milky Way 


29lti«i 


0+62.9 


IC 10 


114t^i 


n+160.5 
'J-0 


LMC 


325+^7^ 


0+56.3 


SMC 


314+167 


0+58.3 



1-2760 
-2740 
■1+4660 
-4610 
-1+2350 
-2430 
-11400 
-11700 
hi 800 
-1780 
-,+302 
-299 
hl07 
-108 
5+327 
-323 

526+264 



5460: 

8820l 

4560J 

22500! 

3400:* 

570j 

2071 

618^ 



7.3 

3.3 

1.1 

3.6 

4.9 









+6.8 

5.1 
+3.7 

2.5 
+5.2 
-1.1 
+2.6 

2.3 
+8.2 

4.9 
32.1 


88.4 


29.6 




-264 







+34.8 



with apparent luminosities > 10^® erg s^^ was calculated on 
the basis of the Hasinger et al. (1998) log(N)-log(S) rela- 
tion, using the distance D and other quantities listed in Ap- 
pendix B (e.g. R25, which was taken from the RC3 catalogue 
of de Vaucouleurs et al. 1991). 

Finally, the number Nulx of ULXs in a galaxy after 
the subtraction of contaminating sources is simply Nulx = 
Nulx, raw — Q- 



2.3 Uncertainties and fitting procedures 

In the following, we consider various sources of error for the 
different considered quantities (SFR, Z, Nulx)- 

For the SFR, we assume that the relative uncertainty 
of each value is ±0.5. This number comes from the analysis 
of the distribution of the measurements reported in Grimm 
et al. (2003): the SFR measurements reported in columns 7, 
8, 9, and 10 of their table 3 were normalized to the 'adopted 
SFR' reported in column 11 of the same table; the distribu- 
tion of the resulting values is bell-shaped and centred around 
1; furthermore, about 68 per cent of them lies between 0.5 
and 1.5. 

We assume that all the metallicity measurements are 
uncertain by 0.1 dex, which is the typical error associated 
with metallicity measurements in HIT regions (see e.g. POl; 
Kennicutt et al. 2003; PT05)"*. 

We note that the uncertainty on the SFR affects the 
error on Nbh in a linear way; instead, because of the compli- 
cated dependence of Nbh upon Z (Nbh changes very rapidly 
with Z for Z ~ 0.3 — OAZq, whereas it is almost indepen- 
dent of Z for Z > 0.5^0 or Z < 0.2^©), the uncertainty 
in Z can result in both very small and very large errors on 
Nbh (this explains the upper limits that can be seen in the 
left-hand panel of Fig. 2). 

In the case of Nulx, raw we adopt Poisson uncertainties, 
as they are affected by small-number fluctuations (Grimm 
et al. 2003). These are estimated through the treatment 
described in Gehrels (1986). Instead, in the case of Q we 
simply assume uncertainties coming from errors upon the 
log(N) - log(S) of Hasinger et al. (1998). 

Finally, the uncertainties on Nulx derive from the com- 
bination of the two above uncertainties. 



■^ The actual error might be larger for several reasons, such as 
uncertainties in the metallicity calibrations, or metallicity fluctu- 
ations within a single galaxy (even at a given radius). 



The analysis presented in the following sections is based 
upon correlation coefficients and power-law fits. Fits were 
performed by minimizing x^ ■ However, since, in most cases, 
the quantities on both axes have uncertainties of similar (rel- 
ative) magnitude, we keep both errors into account through 
the simple procedures summarized in D'Agostini (2005). 

We note that this fact, along with the sizable error bars 
of the variables, tends to produce values of x^ that are 
significantly lower than the number of degrees of freedom 
(hereafter, dof). However, in the low- number regime (that 
we are considering here) the difference between Poissonian 
and Gaussian statistics become significant and hence, when- 
ever possible, we checked our results adopting additional fits 
based on Cash statistics (Cash 1979). In such cases, we fur- 
ther note that i) fits with Cash statistics do not consider 
the errors on the quantity on the a;-axis, and ii) in the cases 
where the variable on the j/-axis is Nulx, we actually per- 
formed a fit to Nulx, raw with the sum of the contamination 
Q to the results of the fitting formula, as required by the 
Cash formalism. The results of the fits based on Cash statis- 
tics are in reasonably good agreement with those of the fits 
based on x^- The main difference is that the Cash statistics 
appears to better constrain the parameters and give smaller 
errors, although this follows in part from neglecting uncer- 
tainties on the a;-axis. In any case, this result reinforces the 
validity of our analysis. 



3 RESULTS 

3.1 Expected number of massive BHs versus 
observed number of ULXs 

Table 2 shows the values of Nbh and ebh (the role of ebh 
is discussed in Appendix C) that we derive for the galaxies 
listed in Table 1, using the models by P98 and by BIO. 
In paper I, we showed that, according to theoretical models 
(Patruno et al. 2005; Blecha et al. 2006), ebh should be of the 
order of 10^^ for the Kroupa IMF, which is consistent with 
most values in Table 2 (especially for the models adopting 
BIO and for galaxies with a large number of ULXs). In Fig. 2, 
we show the observed number of ULXs per galaxy, Nulx, 
as a function of the expected number of massive BHs per 
galaxy, Nbh- 

From Fig. 2 it is immediately evident that there is a cor- 
relation between Nbh and Nulx, when we adopt the models 
from both P98 and BIO. Such correlation can be expressed 
as: 
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Figure 2. Left-hand panel: Number of observed ULXs per galaxy Nulx versus the number of expected massive BHs per galaxy NbHi 
derived using the models from P98. The adopted IMF is a Kroupa IMF (see Appendix A for a comparison with the Salpeter IMF). The 
filled circles are the galaxies listed in Table 1. The solid line is the power-law fit for the entire sample. The dashed line (red on the web) 
is the power-law fit obtained assuming that the index of the power law is /3 = 1. The error bars on both the x— and the y— axis are 1 — ct 
errors. The error bars on Nbh account for the uncertainty on the SFR and on the metallicity (see Section 2.3 for details). Right-hand 
panel: the same as the left-hand panel, but the number of expected massive BHs per galaxy, NbHi has been derived using the models 
from BIO. 
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Figure 3. Nulx versus the SFR. Filled black circles: galaxies 
with metallicity < 0.2 Zq; open circles (red on the web): galaxies 
with metallicity > 0.2 Zq. Solid line: power-law fit for the entire 
sample; dashed line (red on the web): power-law fit obtained as- 
suming that the index of the power law is &i = 1. The error bars 
on both the x— and the y— axis are 1 — tr errors (see Section 2.3 
for details). 



Figure 4. Nulx versus Z. Filled black circles: entire sample; 
solid line: power-law fit. The error bars on both the x— and the 
y— axis arc 1 — a errors (see Section 2.3 for details). 



^ a correlation coefficient with an absolute value 
if the two variables were uncorrelated. 



Pr to obtain 
larger than |r 

The correlation coefficient r is very high, and the prob- 
ability of finding a value larger than |r|, if the two variables 



N 



ULX = 10^NgH, 



(5) 



where /? = 0.80l°;i^ (/3 = 0.851°;^^^ 
(7 = -2.76^" 



-0.12 VM — u-u^-o.is; ^rid 7 — — 2.36_o;g2 
"^0 Ts)' wlifin adopting the models from P98 
(BIO). Table 3 reports the fits, the results of the x^ analysis, 
the Pearson correlation coefficients (r), and the probability 



^ Such probability is based on the assumption that the two vari- 
ables whose correlation is being tested are normally distributed. 
This is likely not the case for our data; however, the value of pr 
can still be retained as an indicator of the strength (or weakness) 
of a correlation. 
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were uncorrelated and normally distributed, is almost zero. 
The values of the x^ associated with the best fits of the 
Nbh — NuLX relation are x ~ 12, with 62 degrees of free- 
dom (dof). The x^ is almost as low when an index /3 = 1.00 
is assumed. According to the F— test (see paragraph 10.2 
in Bevington 1969), the above fit with two parameters (see 
lines 1 and 3 of Table 3) is better than a fit with fixed in- 
dex /3 = 1 (see lines 2 and 4 of Table 3) only at 91 and 
80 per cent confidence level in the case of P98 and BIO, re- 
spectively. This implies that, especially for BIO, the fit with 
two parameters is not a significant improvement. This sug- 
gests that the Nbh — Nulx correlation is a linear relation, 
although deviations from a linear behaviour are possible^. 

Furthermore, the slopes are very similar for the two 
considered models of stellar evolution. The main difference 
between them is the normalization, as the values of Nbh 
obtained assuming P98 are generally a factor of 1.5 — 2 lower 
than the values obtained using BIO. 

Finally, from Tables 2 — 3 and from Fig. 2, it appears 
that the galaxies in our sample that do not host any ULXs 
have very low values of Nbh: according to our model, they 
are consistent with having one or zero ULXs. 

In conclusion. Fig. 2 and Tables 2 and 3 indicate that 
our model predicts a number of massive BHs that correlates 
with the number of ULXs per galaxy, independently of the 
adopted stellar evolution scenario. Such correlation supports 
the hypothesis that (all or a large fraction of) the ULXs are 
connected with massive (rriBH > 25 M©) BHs. 



3.2 Dependence on the SFR and on the 
metallicity 

Our model is based on two key parameters: SFR and metal- 
licity. In this Subsection we discuss their role in more detail. 

3.2.1 Star formation rate 

As mentioned in the Introduction, the existence of a corre- 
lation between SFR and number of bright X-ray sources per 
galaxy is well established (Grimm et al. 2003; Ranalli et al. 
2003; Gilfanov et al. 2004a,b,c). From Fig. 3, it is evident 
that our sample of galaxies follows the same correlation, 
which can be fit as: 



Nt 



10^ 



SFR 

Mo yr- 



(6) 



where S = 0.91 _o it and C = O-l^-oii (with x ~ 18 for 62 
dof, see Table 3), if we consider all the galaxies in the sample. 
We note that the index of the correlation (S) is smaller than 
1, but consistent with 1. Furthermore, the value of the x 
is very low (x^ ~ 18 with 63 dof), even when we assume 
5 — 1.00. According to the F— test, the above fit with two 
parameters (see line 5 of Table 3) is better than a fit with 
fixed index (5=1 (see line 6 of Table 3) only at 48 per cent 



^ The scaling between Nbh and Nulx niight be non-linear be- 
cause Nulx may depend also on properties that do not affect 
Nbh I such as the probability of finding suflBciently massive com- 
panion stars. Thus, a better understanding of the properties of 
HMXBs that can power ULXs is needed, in order to refine our 
model. 



confidence level. This means that the fit with two parameters 
is not an improvement. 

If we impose 5 = 1.00, we obtain Nulx — 
1.20j;jj:i^ SFR. Likewise, Grimm et al. (2003) find an almost 
linear relation between the SFR and the number of X-ray 
sources with Lx > 2 x 10^* erg s~^. Our best-fitting slope 
in the linear relation (1.201q ;[g) is a factor of 2.4±0.4 lower 
than the value (2.9 ±0.23) reported by Grimm et al. (2003); 
given the different luminosity ranges (Lx > 2 x 10^* erg s~^ 
and Lx > lO"'^ erg s~^ in Grimm et al. 2003 and in this 
Paper, respectively), this is expected. In fact, equation (7) 
of Grimm et al. (2003) predicts a factor 3.0 ± 0.5 difference 
between the normalizations of the correlations in the two lu- 
minosity ranges, consistent with our result (the uncertainty 
is due to the error on the slope of the luminosity function 
given in Equation 6 of Grimm et al. 2003). 

We remind also that, in our model, we imposed that 
Nbh scales linearly with the SFR (see equation 3), in 
agreement with the observational behaviour of Nulx. Lines 
13—14 of Table 3 indicate that our fitting procedure infers a 
slope which is consistent with the theoretical value within a 
±0.06 uncertainty. The normalization for the model by P98 
is a factor of ^^^ 1.5 — 2 lower than that of the model by BIO. 

3.2.2 Metallicity 

In Paper I we proposed that there may be an anticorrela- 
tion between the observed number of ULXs (Nulx) and the 
metallicity (Z), but our sample was too small to confirm 
such anticorrelation. Our current sample is more than three 
times larger than the one in Paper I. 

In the current paper, a first argument in favour of the 
role played by the metallicity comes from Fig. 3 (i.e. Nulx 
versus SFR): it is evident that the galaxies with Z > 0.2 Zq 
(open circles) lie mostly below the global fit, whereas the 
galaxies with Z < 0.2 Zq (filled circles) lie preferentially 
above the global fit. This means that a metal-poor {Z < 
0.2 Zq) galaxy tends to have more ULXs with respect to 
a relatively metal-rich {Z > 0.2 Zq) galaxy with the same 
SFR. In quantitative terms, Table 3 shows that, if we split 
our sample according to metallicity and fit the SFR— Nulx 
relation, we obtain similar slopes, but quite different values 
of the normalization for the high-Z sample ((^i = 0.051q'2o) 
and for the low-Z sample (Ci = 0.39to:?3)- The difference 
{0.34:t°oil; or 0.32 ± 0.14 if the slope is fixed to S = 1.00) is 
not consistent with at a significance level slightly below 2 a, 
while it is roughly consistent with the ratio of the expected 
numbers of massive BHs below and above Z ~ 0.2 Zq (see 
Fig. 1). This might be a further indication that the SFR is 
not the only ingredient of the correlation between Nbh and 
Nulx and that the second ingredient might actually be the 
metallicity, as proposed in our model. 

However, we cannot conclude from Fig. 4 and from Ta- 
ble 3 that such anticorrelation exists. There is a very weak 
trend that can be expressed as 



Nulx = lo" {Z/Zq)\ 



(7) 



where rj = -0.21 ± 0.27 and 9 = 0.09 ± 0.20. Such trend is 
not statistically significant, as x^ ~ 86 with 62 dof and the 
correlation coefficient is r ;S — 0.2, with pr = 0.2 (Table 3). 
Similarly to equation (7) , there is no significant anticor- 
relation, in our model, between Nbh and Z (lines 15-16 of 
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Table 3. Parameters of the power-law fits and x 



X— axis 


y~ axis 


Model 


Sample ^ 


Index ''•<= 


Normalization "^ 


xVdof <* 


r (Pr) ^ 


Nbh 


NuLX 


P98 


all 


0.80+° :[^(0.86± 0.07) 


-2.3610-45 (_2.64± 0.28) 


8.7/62 


0.90 (2x10^23) 


Nbh 


NuLX 


P98 


all 


1.00(1.00) 


-3.11 ±0.07(-3.17± 0.04) 


10.0/63 


0.90 (2x10-23) 


Nbh 


NuLX 


BIO 


all 


0.85lO-^^{0.90±0.07) 


-2.76l°'^^(-3.00±0.30) 


11.1/62 


0.93 (1x10-27) 


Nbh 


NuLX 


BIO 


all 


1.00(1.00) 


-3.36 ±0.07(-3.41± 0.04) 


11.8/63 


0.93 (1x10-27) 


SFR 


NuLX 


- 


all 


0.91+O;^^(0.90±0.08) 


0.13l[;;JO(Q09±o.06) 


17.7/62 


0.88 (4x10-22) 


SFR 


NuLX 


- 


all 


1.00(1.00) 


0.08 ±0.06(0.02 ±0.04) 


17.8/63 


0.88 (4x10-22) 


SFR 


NuLX 


- 


lowZ 


0.75+°;^°(0.80±0.08) 


0.39l°;°^(0.32±0.07) 


4.4/22 


0.93 (8x10-11) 


SFR 


NuLX 


- 


highZ 


0.83t°;^^(0.94±0.15) 


0-05-a2o(-0-06±0.10) 


6.4/38 


0.88 (6x10-") 


Z/Zq 


NuLX 


- 


all 


-0.21 ±0.27(-0.55± 0.15) 


0.09 ±0.20(-0.04± 0.12) 


86.0/62 


-0.16 (2x10-1) 


Z/Zq 


NuLX 


- 


all 


0.00(0.00) 


0.23 ±0.05(0.36 ±0.04) 


86.6/63 


-0.16 (2x10-1) 


Z/Zq 


NuLx/SFR 


- 


all 


-0.55 ± 0.23 


-0.37 ±0.18 


10.4/62 


-0.30 (2x10-2) 


Z/Zq 


NuLx/SFR 


- 


all 


0.00 


-0.03 ±0.07 


14.7/63 


-0.30 (2x10-2) 


SFR 


Nbh'' 


P98 


all 


0.96 ±0.06 


3.19 ±0.04 


13.8/62 


0.82 (7x10-1^) 


SFR 


Nbh' 


BIO 


all 


0.97 ±0.05 


3.44 ± 0.04 


6.3/62 


0.95 (6x10-33) 


Z/Zq 


Nbh' 


P98 


all 


-0.19 ±0.29 


1 41+0.23 


153.9/62 


-0.23 (7x10-2) 


Z/Zq 


Nbh' 


BIO 


all 


0-05l«:f, 


1.85 ±0.24 


183.2/62 


-0.11 (4x10-1) 


Z/Zq 


Nbh/SFR f 


P98 


all 


-0.60 ±0.07 


2.79 ± 0.05 


9.4/62 


-0.96 (2x10-37) 


Z/Zq 


Nbh/SFR f 


BIO 


all 


n 34+0-02 


3.22 ± 0.04 


17.0/62 


-0.98 (3x10-18) 



The SFRs used by the fitting procedure are in units of Mq yr" . ^ The sample adopted for the fits is referred to as 'all' when all the 
galaxies in Table 1 are considered, and as lowZ (highZ) when only the galaxies with Z < 0.2 Zq (Z > 0.2 Zq) are considered. ^ When 
the index is equal to 1.00 or to 0.00, without error, it means that it has been fixed to such value. '^ Values in parenthesis are the fitting 
parameters obtained with Cash statistics (where applicable), rather than x2. such results do not take into account the uncertainty on 
the quantity on the x-axis. "^ 'x^ /doi is the x^ divided by the degrees of freedom (dof). "^ r is the Pearson correlation coefficient, and Pr 
is the probability of finding a value larger than jr|, if the two variables were uncorrelated and normally distributed. The fits listed in 
rows 13-18 refer entirely to the theoretical model. We report them only for a comparison between the statistical fluctuations of the 

theoretical calculations and those of the observational data. 
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Figure 5. From left to right: Nulx/SFR versus Z, Nbh/SFR versus Z for P98, Nbh/SFR versus Z for BIO. Filled black circles: entire 
sample. Solid lines: power-law fit. The error bars on both the x— and the y— axis are 1 — a errors. Central and right-hand panel: the 
error bars on Nbh/SFR account for the uncertainty on the metallicity (see Section 2.3 for details). 



Table 3) . In particular, the values of the x^ for the best fits 
are 154 and 183 for P98 and BIO, respectively, with 62 dof 
(see lines 15 and 16 of Table 3). The correlation coefficients 
are low: r= —0.2 and = —0.1 for P98 and BIO, respectively. 
This happens because, in our model, the dependence of Nbh 



on the SFR dominates with respect to the dependence of 
Nbh on the metallicity. 



To better highlight the possible effects of metallicity, 
we plot the number of ULXs per galaxy Nulx normalized 
to the SFR versus the metallicity (left-hand panel of Fig. 5). 
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The fit of tlie left-fiand panel in Fig. 5 gives tfie following 
results. 



NuLx^^|^')=10-(Z/Zer, 



where ti = -0.55±0.23 and ki = -0.37±0.18, with x^ 



(8) 
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for 62 dof (and r= —0.3). Thus, there is marginal evidence of 
an anticorrelation between Z and Nulx/SFR. According to 
the F— test, the above fit with two parameters (see line 11 of 
Table 3) is significantly better (at a 96 per cent significance 
level) than a fit with fixed index ti = (i.e. the case of no 
correlation) and with only one free parameter (ki = — 0.03± 
0.07, line 12 of Table 3). This result supports the hypothesis 
of an anti-correlation between Nulx/SFR and Z. 

As in the case of the comparison with Nulx, the role 
of metallicity in our model can be better appreciated if we 
remove the contribution from the SFR. In fact, the plot 
of Nbh/SFR versus Z (central and right-hand panels of 
Fig. 5) shows quite clearly the effect of metallicity in equa- 
tion (1). We know from equation (1) that the relation be- 
tween Nbh/SFR and Z is not a power-law. However, in order 
to compare it with equation (8), we can approximate it as: 



Ne 



Mq yr 
SFR 



10"^ {z/z^y^ 



(9) 



where t2 = -0.60 ± 0.07 i~0-S4,to'ol) and Ka = 2.79 ± 0.05 
(3.22 ± 0.04) for P98 (BIO). The statistical significance of 
such correlations is given in lines 17—18 of Table 3. ti is 
consistent with t2 in the case of P98 and marginally consis- 
tent with i2 in the case of BIO. Thus, we can conclude that 
the anticorrelation in the central and right-hand panels of 
Fig. 5 is consistent with that in the left-hand panel of the 
same Figure, although the error bars are quite large. 

The fact that the anticorrelation between Z and Nulx 
only emerges after removing the effect of the SFR (i.e. con- 
sidering NuLx/SFR versus Z) clearly shows that the SFR 
(and not the metallicity) is the dominant factor in determin- 
ing the number of ULXs in a given galaxy^. However, the 
metallicity plays a crucial role in enabling the formation of 
ULXs, as massive BHs can be produced only in star-forming 
regions where Z is sufficiently low. 

We emphasize that, in our model, the effect of the 
metallicity on the formation of massive BHs is mainly a 
threshold effect: in the models by BIO (see their fig. 1 and 
their equation 11) massive (ttibh > 25 Mq) BHs form only 
for Z < 0.47 Zq. If the metallicity Z is below this threshold, 
there is quite a small spread in rriprog (defined in equation 1) 
and thus in Nbh for various metallicities (at fixed SFR): 
adopting the model by BIO, rriprog ranges from 68 Mq for 
Z = 0.04^0 (the Antennae) to 96 M© for Z = 0.45^0 
(NGC 2442). Furthermore, as Nbh scales linearly with the 
SFR, the effect of metallicity on Nbh is small with respect 
to the effect of the SFR, for Z < 0.47 Zq. 

Instead, if Z > 0.47 Zq, massive BHs cannot form from 
the direct collapse of massive stars. This might be a prob- 
lem for our model, as ULXs exist also in a few galaxies 
with Z > 0.47 Zq: NGC 4501 (see Table 1) is the only 
galaxy, for which we found SFR, Z and X-ray data, that has 



^ The importance of the SFR is somewiiat amplified by the larger 
spread in the SFR values with respect to the metallicity ones. 



Nulx, raw > and Z > 0.47 Zq . In such a case, the contami- 
nation from foreground/background sources (Q — 0.43lo Jg) 
is not completely negligible; nonetheless it is insufficient to 
explain the 2 observed ULXs. However, there are several 
reasons for which this might happen. 

First, as the metallicity is uncertain by ~ 0.1 dex, val- 
ues of Z below 0.47 Zq are well within 1 a error ranges for 
NGC 4501 (1- a lower limit on Z is 0.39^0 for NGC 4501). 
Thus, the upper limits of Nbh for NGC 4501 are larger than 
zero. 

Second, in the case of large spiral galaxies (such 
as NGC 4501) our choice of taking the metallicity at 
R — 0.73 i?25, although justified, is rather simplistic: in 
NGC 4501 the metallicity gradient is such that Z goes below 
0.4 Zq ai R < R25, so that there exists a small region where 
massive BHs might form, even in the P98 model. Thus, 
uncertainties in the metallicity determination and spatial 
abundance fluctuations might further help in this respect. 

Third, in this Paper we considered only BHs with 
m-BH > 25 M0 as engines of the ULXs, but slightly less 
massive BHs might still power many low- luminosity ULXs, 
provided that there is a certain level of beaming or super- 
Eddington accretion (e.g. King 2008): for example, if the 
metallicity is Z ^ 0.57 Zq BHs with iubh ~ 20 M© can still 
form (BIO). It is worth mentioning that NGC 4501 hosts rel- 
atively faint ULXs, that can be explained without invoking 
massive BHs, nor strong beaming/super-Eddington accre- 
tion. 

Finally, the metallicity needed in our model is that of 
the molecular clouds before the pollution from the first su- 
pernovae associated with the parent cluster, as massive stars 
(> 40 M0) collapse into BHs before the explosion of such 
supernovae. Thus, the metallicity measured today is likely 
higher than the value we should consider in our model, es- 
pecially for post-starburst galaxies. 



4 CONCLUSIONS 

Low-metallicity {Z ^ 0.4 Zq) massive ( ^ 40 M0) stars are 
expected to produce massive remnants (25 < thbh/M© < 
80, H03, BIO) at the end of their evolution. Such massive 
BHs might power a large fraction of the observed ULXs 
in low-metallicity galaxies. In this Paper, we derived the 
number of massive BHs (Nbh) that are expected to form in 
a galaxy, via this mechanism, in the same time period that 
they could have a massive ( ;J, 15 Mq) donor companion. 
We find that Nbh correlates well with the observed number 
of ULXs per galaxy (Nulx). The slope of such correlation 
does not depend significantly either on the IMF or on the 
adopted stellar evolution model. The IMF and the stellar- 
evolution models affect only the normalization of Nbh (the 
spread is generally less than a factor of 2). We stress that 
the stellar-evolution models adopted in this Paper neglect 
some important effects, such as the rotation and the possible 
influence of binary evolution. The final mass of the remnant 
is likely affected by the fact that the massive progenitor of 
the BH is fast rotating or that it resides in a binary, where 
additional mass loss is possible (see e.g. H03). Accounting 
for the probability that the progenitor of the massive BH is 
in a binary system likely introduces additional uncertainties. 
In addition, the model considered here does not in- 
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elude the possibility of pair instability supernovae (PISNs) . 
PISNs are predicted to occur at very low metallicity, in the 
case of very massive stars (> 140 M©), and lead to the 
complete disruption of the star (HW02; H03). Pulsational 
pair instability may occur also for smaller stellar masses 
(100 — 140 M0) and/or for larger metallicity, but it does 
not lead to a PISN, and a massive BH can form via direct 
collapse. Thus, PISNs probably do not play a role for stars 
with metallicity Z > 0.01 Zq and mass < 140 M© (H03). 
On the other hand, even assuming (as a strong upper limit) 
that all stars with mass > 100 Mq do not leave any rem- 
nant, due to a PISN, our estimates of Nbh change by less 
than 10 per cent. 

The model described in this Paper is consistent with the 
observed correlation between the number of ULXs and the 
SFR, as well as with the fact that ULXs are found prefer- 
entially in low-metallicity environments. Furthermore, this 
model is a natural extension of the scenario described by 
Grimm et al. (2003). In fact, Grimm et al. (2003) find a cor- 
relation between the SFR and the number of X-ray sources 
(not necessarily ULXs) per galaxy, and they explain it with 
the correlation between the SFR and the number of HMXBs 
powered by stellar BHs. In this Paper we suggest that this 
correlation still holds for the ULXs, as the ULXs (or most 
of them) can be powered by massive BHs formed by the 
collapse of massive metal-poor stars. 

Our model predicts the existence of a dependence of 
NuLX on Z (and the data suggest it, too), when the dom- 
inant effect due to the SFR is removed. Unfortunately, the 
statistical uncertainty of such dependence is still quite high, 
due to the dearth and to the inhomogeneity of the data. 

In particular, there is a large inhomogeneity in the 
metallicity measurements. Moreover, the metallicity needed 
in our model is that of the molecular clouds before the pol- 
lution from the first supernovae associated with the parent 
clusters, as very massive stars collapse into BHs before such 
supernovae. Thus, the metallicity measured today is likely 
higher than the value we should consider in our model. Only 
for some types of galaxies, such as the ring galaxies, where 
the star-formation history has a clear connection with the 
geometry of the system, is it possible to measure a pre- 
starburst value of Z, suitable for our purposes. A possible 
way to reduce this problem and to check our model is to take 
new measurements of the local metallicity in the neighbour- 
hoods of the observed ULXs, or even in the nebula associated 
with the ULXs (Ripamonti et al., in preparation). 
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APPENDIX A: COMPARISON BETWEEN 
KROUPA AND SALPETER IMF 

In this paper we assumed the Kroupa IMF (Kroupa 2001), 
a moderately top-heavy IMF. In this Appendix we show 
how sensitive our model is to the adopted IMF, by using 
the Salpeter IMF (Salpeter 1995) to derive Nbh. Table Al 
and Fig. Al show the results obtained adopting the Salpeter 
IMF. Basically, all the values of Nbh calculated with the 
Salpeter IMF are a factor of 2 lower than the ones ob- 
tained with the Kroupa IMF. However, the slope of the fit 
in equation (5) does not change significantly with a differ- 
ent IMF. Adopting the Salpeter IMF (Fig. Al), for the fits 
in equation (5) we find /3 = O.SOlg:!! {l3 = Q.SSto.lt) and 



7: 



-2.12l»;«(7^ 



-2.50lg;7Q), corresponding to x^ = 8.7 



(X = 11-1) with 62 dof, for the models by P98 (BIO). As for 
the Kroupa IMF, fixing /3 = 1, we still obtain good values 
of X^ (X^ = 10-0 and 11.7 with 63 dof, for P98 and BIO, re- 
spectively). Thus, only the normalization of Nbh is affected 
by a different choice of the IMF. 



APPENDIX B: DATA COLLECTION AND 
ELABORATION 

In this section we provide details about how we collected and 
calculated the values reported in Table 1, that represent the 
observational basis of this paper. 



Bl Star Formation Rate 

As we already mentioned, the SFR reported in Tables 1 
comes from either ultra-violet (UV), Ha, far-infrared (FIR), 
or radio measurements. If more then one independent mea- 
surement is available, we generally average the correspond- 
ing SFRs, keeping into account the properties of the galaxy, 
of the observational data, and of calibrations. 



In particular, for galaxies included in the UV cata- 
logue by Mufioz-Mateos et al. (2007), we derive SFR= 
27rr^ EsFRO [1 - (H- not/ra) exp (-rtot/ra)], where Esfrq 
is the central SFR surface density, r^ is a characteristic 
length-scale and rtot is the external radius of the galaxy 
(from the integration of equation Bla in Mufioz-Mateos et 
al. 2007). When the total Ha luminosities are available, we 
adopt the correlation in Kennicutt (1998; hereafter K98), 
i.e. SFR= L(Ha)/1.26 x 10^^ erg s'^ [Mq yr'^]. Similarly, 
when the FIR luminosity is available, we apply the equation 
SFR= L(FIR)/2.2 x 10^^ erg s"^ [M© yr'^] in K98. When 
a radio measurement (at 1.4 GHz) is available, we then use 
equation (6) of Bell (2003). Finally, if a galaxy belongs to 
the sample by Grimm et al. (2003), we generally use their 
adopted values. 



B2 Metallicity 

The existence of a metallicity measurement is the most re- 
strictive among the criteria for inclusion in our sample, as 
metallicities are often unavailable. 

For most galaxies in the sample we use measurements of 
oxygen abundance (based on oxygen line intensities from HII 
regions) as proxies for metallicity. Instead, when the spec- 
tra of HII regions are unavailable, we use X-ray metallicity 
estimates (from papers on X-ray measurements of ULXs), 
although they are much less accurate. In the rest of this 
subsection we provide extra details about the measurements 
from HII regions. 

B2. 1 Conversion of line intensities to metallicities 

When line intensity measurements include the weak 
A=4363A OIII line (e.g. in the case of HII regions of the 
Cartwheel galaxy from Fosbury & Hawarden 1979), we sim- 
ply use the oxygen abundances reported in the literature: 
this is possible because the measurement of the 4363A line 
greatly simplifies the conversion of line intensities into abun- 
dances, and the conversion procedure did not change much 
since it was first established (see e.g. Pagel et al. 1992; and 
the discussion in PT05). Instead, when only the most in- 
tense OH (A=3727A) and OIII (A=4959,5007A) lines are 
observed, the oxygen abundances in different papers are of- 
ten based on different calibrations (see the discussion in POl, 
and in PT05), that lead to significant offsets in the results. 
For this reason, we went back to the observational data and 
applied the PT05 calibration whenever this was possible. We 
made an exception for galaxies where the metallicity value in 
the literature had been calculated with the POl calibration 
(such as the ones in the compilation by Pilyugin, Vilchez 
& Contini 2004; hereafter PVC04); this is because the two 
calibrations are quite close to each other (the PT05 calibra- 
tion is essentially a revision of the POl calibration), and the 
difference (typically amounting to ~ 0.05 dex for single HII 
regions; see the top panel of Fig. 7 from PT05) should not 
have a large effect. 

We note that both the PT05, and the POl calibrations 
generally provide significantly lower abundances with re- 
spect to previous calibrations (e.g. Edmunds & Pagel 1984), 
but are currently considered the most accurate (see e.g. Ken- 
nicutt, Bresolin & Garnett 2003, who actually suggest that 
abundances might be even lower). 
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Table Al. Values of Nbh and ebh for the galaxies listed in Table 1, assuming the models by P98 (first and second column) and by BIO 
(third and fourth column). We adopt a Salpeter IMF. 



Galaxy 



Nbh (P98) ebh/IO-* (P98) Nbh (BIO) eBn/lQ-^ (BIO) 
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Figure Al. Number of expected massive BHs per galaxy, derived using the models from P98 (left-hand panel) and from BIO (right-hand 
panel), versus the number of observed ULXs per galaxy. The adopted IMF is a Salpeter IMF. The filled circles are the galaxies listed in 
Table 1. The error bars on both the x— and the y— axis are 1 — a errors. The error bars on Nbh account for the uncertainty on the SFR 
and on the metallicity (see Section 2.3 for details). The solid lines are the power-law fits, for the entire sample. The dashed lines (red on 
the web) are the power-law fits. 



Table Al - continued. 



Galaxy 



Nbh (P98) ebh/IO-* (P98) Nbh (BIO) eBn/lO-* (BIO) 



Circinus 
NGC 6946 
PGC 68618 
NGC 7714/7715 
NGC 7742 
Milky Way 
IC 10 
LMC 
SMC 



f920 
-920 
2270+1250 
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142 
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+79 
142 
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2.3 
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lO.O 
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When the metallicity gradient is available for a given 
galaxy, we list in Table 1 and adopt in our calculations the 
value of the metallicity at r = 0.73 i?25, where i?25 is the 
isophotal or Holmberg radius. In fact, r = 0.73 7?25 is the 
average distance of the observed ULXs from the centre of 
the host galaxy, in a sample of spiral galaxies (fig. 15 in Liu 
et al. 2006; the averaging process took contamination into 
account). 

The oxygen abundance is usually derived in terms of 
12-|-log(0/H); this can be easily converted to the units used 
in Table 1 (where the metallicity Z is expressed as a fraction 
of solar metallicity Zq), by assuming that 12-|-log(0/H)Q = 
8.92, which corresponds to Zq — 0.02 (the value commonly 
used in stellar evolution studies). 



B3 Distance 

For the 52 galaxies in the LB05 catalogue we use the dis- 
tance reported there. For the MW we simply indicate the 
distance to the Galactic centre (8.5 kpc). For the remain- 
ing 11 galaxies we use the distances reported by the NASA 



Extragalactic Database (NED). In the 9 cases where it is 
possible, we take the mean redshift-independent distance; 
for the remaining 2 galaxies we derive distances from the 
redshift, assuming Ho = 73kms^^ Mpc"'^. 



B4 Contaminating sources 

B4-1 Further required data 

Estimating the number of contaminating sources requires 
the use of several quantities that are not reported in Table 1: 
the galaxy angular size (major and minor axis R25 and r25 
of the 25 magarcsec"^ isophote), the sky area ^obs actually 
covered by the X-ray observations, and the total galactic HI 
column density (Nh). 

For galaxies in the LB05 catalogue, we use the R25, 
r25 and A^n values reported there. For the others, we adopt 
the references quoted in LB05: the RC3 catalogue (de Vau- 
couleurs et al. 1991) for angular sizes, and the paper of 
Dickey & Lockman (1990; hereafter DL90) for the weighted 
average of Ah. We make an exception for the MW and the 
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Magellanic Clouds, where such quantities do not really mat- 
ter, since the number density of the expected contaminating 
sources is extremely low, and there are no ULXs inside such 
galaxies. 

The sky area j4obs was taken from the papers report- 
ing the X-ray observations. When no specific information is 
provided, we assume an area of 314 arcmin^ if the data were 
taken with ROSAT (this corresponds to the area within 10' 
from the instrument's axis, where the sensitivity is reason- 
ably constant and close to the maximum value; see LB05), or 
an area of 70.6 arcmin^ (i.e. the entire 8.4' x 8.4' field of view 
of the S3 detector) if the data were taken with Chandra. 

B4-2 Contamination estimate 

In order to estimate the number of background or foreground 
contaminating sources, we followed a procedure very similar 
to the one reported by Liu et al. (2006). 

Since the bulk of the galaxies in our sample was ob- 
served with ROSAT, we used the log(N)-log(S) reported in 
Hasinger et al. (1998). We converted the minimum lumi- 
nosity of an ULX, Lum ~ lO'^^ergs^^, into the apparent 
flux limit Fiini = Liini/(47rD^); then, we converted Fum into 
another flux 5iim that can be used in the Hasinger et al. 
(1998) log(N)-log(S) relation. Such conversion takes into ac- 
count (i) the different assumptions on the observed band 
(the Hasinger et al. 1998 flux refers to the 0.5-2.4 keV band, 
whereas most references for ULXs provide luminosities in 
other bands - e.g. 0.3-8.0 keV band); (ii) the different as- 
sumption over spectral slopes (Hasinger et al. 1998 assumes 
that all sources have photon index 2; whereas here we as- 
sume a photon index 1.7, as done by LB05); (iii) the absorp- 
tion from the galactic A'^e. In most cases, such conversion 
amounts to a reduction by a factor of 2-3 (the exact value 
depending on the band used for the X-ray observations, and 
to a lesser extent on the galactic Ah for each specific galaxy) 
of -Film- Finally, we integrate the Hasinger et al. (1998) dif- 
ferential log(N)-log(S), and find the expected surface num- 
ber density q of contaminating sources (i.e. of sources with 
apparent flux larger than 5iim). 

In order to determine the expected contamination in 
a specific galaxy, we combined q with the size of the 
galaxy (and of the observed area), taking into account 
also the radial distribution of ULXs. The exact proce- 
dure is the following. First of all, we calculate the ob- 
served area with deprojected distance < R25 from the 
centre Ai = min(Aobs,7r-R25»'25), and the observed area 
with deprojected distance between R25 and 2R25, A2 = 
min(^obs — Al, 37r_R25''25). This assumes that the field of 
the X-ray observation was centred on the host galaxy, and 
that such field is circular. Then, we obtain the number of ex- 
pected contaminating sources in the two regions, Qi — Aiq 
and Q2 = A2q. These must be compared to the numbers Ni 
and N2 of observed ULXs in the corresponding areas (obvi- 
ously, Ni -I- N2 = NuLX,raw), siucc the number of contami- 
nating sources cannot be larger than the number of observed 
sources. Then, we estimate the actual number of contami- 
nating sources to be 



i(Qi,Ni)-hmin(Q2,N2 



(Bl) 



This expression automatically sets the number of contami- 
nating sources at for galaxies where no ULX was observed. 



For the 10 galaxies where the original X-ray obser- 
vations provide a value of the number of ULXs that was 
already cleared from contaminating sources (e.g. because 
one or more of the candidate ULXs in a galaxy have been 
identified as background objects), we simply assume that 
NuLx.raw — NuLx, i.e. that Q = 0. 

The total number of expected contaminating sources 
in the 54 galaxies where we apply our estimates is 14.60, 
whereas the sum of NuLx.raw in such galaxies is 112; the 
contamination fraction is then ~0.13, which is smaller than 
those found in previous papers (25 — 44 per cent, Swartz 
et al. 2004; Ptak & Colbert 2004; Liu et al. 2006; Lopez- 
Corredoira & Gutierrez 2006). There are two reasons for 
this discrepancy. First, we do not consider elliptical galax- 
ies, where the contamination is stronger: for example, Liu 
et al. (2006) estimate a contamination fraction ~ 1 for early 
type galaxies, whereas for late-type galaxy such fraction is 
only ~ 0.2. Second, our sample includes a smaller fraction of 
distant galax;ies than most ULXs catalogues (for example, 
in the LB05 catalogue the fraction of galaxies with distances 
> 20Mpc is ~ 0.4, whereas in our sample it is ^ 0.1): this 
is relevant because distant galaxies tend to have higher con- 
tamination levels than nearby galaxies*. 



B5 A possible contamination from Population II 
ULXs 

ULXs in elliptical galaxies might be connected to the minor 
(but non-negligible) fraction of ULXs that is claimed (e.g. 
Colbert et al. 2004; Brassington et al. 2005) to be associated 
with old stellar populations (Population II ULXs). Since old 
populations are present also in spiral galaxies. Population II 
ULXs might be present also in our sample. It is quite difficult 
to estimate their number. 

Liu et al. (2006) estimate that the number of ULXs 
in early-type and elliptical galaxies can be completely ex- 
plained by contaminating sources (i.e. without requiring 
Population II ULXs): if so, the number of ULXs associated 
to old stellar populations should be very low (or even neg- 
ligible) also in spiral/irregular galaxies. On the other hand, 
Colbert et al. (2004) suggest that Population II ULXs repre- 
sent a fraction ~ 0.2 of ULXs in spirals. However, the results 
by Colbert et al. (2004) and by Liu et al. (2006) are not nec- 
essarily in conflict. In fact, l-cr upper limits for early-type 
galaxies® from Liu et al. (2006) are not far from the Col- 
bert et al. (2004) results. Since the fraction of population II 



° The distance D affects of the contamination Q through the 
angular surface of a galaxy, and the corresponding limit flux 5ii m . 
Both scale as D~^. However, the Hasinger et al. (1998) log(N) — 
log(S) implies that for D ^ 30 Mpc the decrease in S'lim produces 
a steep increase in the surface density of contaminating sources 
{q oc D^*), and Q oc D^*; instead, for sources at large distances 
the two dependencies almost cancel out, and Q oc D"'''^. 
^ In their Section 3.2 and Tables 1 and 2, Liu et al. (2006) esti- 
mate the occurrence frequencies of ULXs. They find 0.72 it 0.11 
ULXs per spiral galaxy, 0.02 ± 0.11 ULXs per early- type galaxy, 
and —0.15 ± 0.13 ULXs per elliptical galaxy. The ratio is then 
0.03 ± 0.15 if wc consider early-type galaxies, and —0.21 ± 0.19 if 
wc consider ellipticals. The ratio for early-types can be taken as 
an upper limit on the fraction of Population II ULXs in spirals, 
since the old population of early- typc/cUiptical galaxies of the Liu 
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ULXs in spiral galaxies is small for any of these estimates, 
we neglect their contribution to our ULX sample. 



B6 References and details for the single galaxies 

Here we give the detailed references for the data used to 
compile Table 1, grouped by galaxy. 

The Cartwheel (ESO 350-G 040): average SFR 
from Appleton & Marston (1997, Ha data) and from Mayya 
et al. (2005, radio measurement); Z from Fosbury & Hawar- 
den (1977), spectra of HIT regions in the outer ring; X- 
ray sources from Wolter & Trinchieri (2004). This galaxy 
is not in the LB05 catalogue. Distance from the NASA 
Extragalactic Database (hereafter NED), assuming Ho = 
73kms"^Mpc^^ and T/(LocalGroup) = 9048 km s"^ Angu- 
lar size and axis ratio from the RC3 catalogue; A'^e from 
DL90. 

NGC 253: SFR from Grimm, Gilfanov & Sunyaev 
(2003); Z from Webster & Smith (1983), that reports the 
measurements of line intensities for OII[A3727 + A3729] and 
OIII[A4959 + A5007], from which we derive the metallicity 
with the formula by POl; X-ray from LM05. 

NGC 300: SFR from Helou ct al. (2004), Ha mea- 
surement, using the correlation by K98; metallicity gradient 
from PVC04; X-ray from LB05. 

NGC 598 (M 33): SFR from ultra-violet (Mufioz- 
Mateos et al. 2007); metallicity gradient from spectra of HIT 
regions (PVC04); X-ray from LM05. This galaxy is not in 
the LB05 catalogue. Distance from NED (mean redshift- 
independent distance). Angular size and axis ratio from the 
RC3 catalogue; A''h from DL90. 

NGC 628 (M 74): SFR from Grimm et al. (2003); 
metallicity gradient from spectra of HIT regions (Pilyugin et 
al. 2002); X-ray from LM05. This galaxy is not in the LB05 
catalogue. Distance from NED (mean redshift-independent 
distance). Angular size and axis ratio from the RC3 cata- 
logue; A^H from DL90. 

NGC 1058: SFR from Ha (Kennicutt& Kent 1983), 
using the correlation by K98; metallicity gradient from 
Moustakas & Kennicutt (2006a), based on the PT05 cali- 
bration; X-ray from LM05. This galaxy is not in the LB05 
catalogue. Distance from NED (mean redshift-independent 
distance). Angular size and axis ratio from the RC3 cata- 
logue; A^H from DL90. 

NGC 1073: SFR from Ha (Martin & Friedh 1997), us- 
ing the correlation by K98; metallicity gradient from spectra 
of hi: regions (Dors & Copetti 2005); X-ray from LB05. 

NGC 1291: SFR from ultra-violet (Muiioz-Mateos et 
al. 2007); Z from X-ray measurements (Irwin, Sarazin &: 
Bregman 2002); X-ray from LM05. 

NGC 1313: SFR from Ha (Ryder & Dopita 1994); Z 
from spectra of HII regions (POl); X-ray from Colbert et al. 
(1995). 

NGC 1365: average SFR from mid- and far-infrared 
measurements (Lonsdale & Helou 1985; Strateva & Komossa 
2009a, 2009b); metallicity gradient from spectra of HII re- 
gions (PVC04); X-ray from Strateva & Komossa (2009a). 



ct al.(2006) sample is typically larger than the old population in 
the spiral/irregular galaxies of the same sample. 



IC 342 (PGC 13826): SFR from Grimm et al. (2003); 
metallicity gradient from spectra of HII regions (PVC04); 
X-ray from LM05. 

NGC 1566: SFR from far-infrared (SFRfir = 
3.24Mq yr~^, Bell & Kennicutt 2001, using the correlation 
by K98) ; Z from spectra of HII regions (Hawley & Phillips 
1980); X-ray from LB05. 

NGC 1705: average SFR from two different Ha mea- 
surements (SFRHa = 0.063 Moyr^^ Gil de Paz, Madore 
& Pevunova 2003; SFRhc = 0.112 M0yr"\ Hunter & 
Elmegreen 2004), using the correlation by K98; Z from Lee 
& Skillman (2004), spectra of HII regions; X-ray from LB05. 

NGC 2366: average SFR from Ha (SFRh,, = 
0.125 Meyr^^ Hunter & Elmegreen 2004, using the cali- 
bration by K98, 8-1000 pm (SFRir = 0.03 Mgyr^^ Hop- 
kins, Schulte-Ladbeck & Drozdovsky 2002, using the cal- 
ibration by K98 and radio (SFRi.4ghz = 0.069 Moyr^^ 
Condon, Cotton & Broderick 2002, using the calibration by 
Bell 2003); Z from spectra of HII regions (PVC04); X-ray 
from LB05. 

NGC 2403: SFR from Gilfanov, Grimm & Sunyaev 
(2004a); metallicity gradient from spectra of HII regions 
(PVC04); X-ray from LM05. 

NGC 2442: average SFR from ultra-violet (SFRuv = 
3.86 Mq yr"\ Muiioz-Mateos et al. 2007) and Ha (SFRh^ = 
5.26MQyr~^, Helmboldt et al 2004, using the calibration 
by K98); metallicity gradient from spectra of HII regions 
(PVC04); X-ray from LB05. 

Holmberg II (Arp 268, PGC 23324): SFR from 
Ha (Walter et al. 2007), adopting the relation in K98; Z 
from spectra of HII regions (Walter et al. 2007); X-ray from 
Dewangan et al. (2004). 

NGC 2903: average SFR from far-infrared (SFRfir = 
1.8MQyr^^, Ranalli et al. 2003, using the calibration by 
K98; SFRfir = 1.5 Mq yr~\ Bell 2003, using the cahbration 
by K98) and radio (SFRi.4GHz = 2.6M0yr"\ Ranalh et 
al. 2003, using the calibration by Bell 2003; SFRi.4ghz = 
1.6 Mq yr~^, Bell 2003); metallicity gradient from spectra of 
HII regions (Pilyugin, Thuan & Vilchez 2006); X-ray from 
LB05. 

NGC 3031 (M 81): SFR from ultra-violet (Mufioz- 
Mateos et al. 2007); metallicity gradient from spectra of HII 
regions (PVC04); X-ray from LM05. 

NGC 3049: SFR from radio (Bell 2003); average Z 
from Gusova, Izotov & Thuan (2000) and Kehrig, Telles & 
Cuisinier (2004); X-ray from LB05. 

PGC 30819 (IC 2574, UGC 05666): average SFR 
from Ha (SFRh^ = 0.076 MQyr-\ Miller & Hodge 1994, 
using the calibration by K98) and radio (SFRi.4ghz = 
0.039 Mq yr~^, Condon et al. 2002, using the calibration by 
Bell 2003); Z from spectra of HII regions (Miller & Hodge 
1996); X-ray from LB05. 

NGC 3310 (Arp 217): SFR from Grimm et al. (2003); 
Z from spectra of HII regions (Pastoriza et al. 1993) , metal- 
licity from spectra of HII regions; X-ray from LM05. 

NGC 3395/3396 (Arp 270): SFR from radio (Con- 
don et al. 2002), using the calibration by Bell (2003); Z from 
integrated spectrum (Kennicutt 1992); X-ray from Brassing- 
ton et al. (2005). 

PGC 35286 (UGC 06456): SFR from Ha (Hunter 
& Elmegreen 2004), using the calibration by K98; Z from 
spectra of HII regions (Pilyugin 2001b); X-ray from LB05. 
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PGC 35684 (UGC 06541, Mkn 178): SFR from Ha 
(Hunter & Elmegreen 2004), using the calibration by K98; 
Z from spectra of HII regions (Pilyugin 2001b); X-ray from 
LB05. 

NGC 3738 (Arp 234): average SFR from Ha 
(SFRhc« = O.O47M0yr"\ Hunter & Elmegreen 2004) and 
radio (SFRi.4GHz = 0.029 M0yr"\ Condon et al. 2002, us- 
ing the calibration by Bell 2003); Z from spectra of HH 
regions (Martin 1997); X-ray from LB05. 

NGC 3972: average SFR from Ha (SFRho = 
0.15 Mq yr~^, Moustakas & Kcnnicutt 2006b, using the cali- 
bration by K98), far-infrared (SFRfir = 0.38 Mq yr"^ Bell 
2003, using the calibration by K98) and radio (SFRi.4ghz = 
O.36M0yr-\ Bell 2003); Z from spectra of HH regions 
(PT05); X-ray from LB05. 

The Antennae (NGC 4038/4039, Arp 244): SFR 
from Grimm et al. (2003); X-ray estimate of Z (Fabbiano 
et al. 2004); X-ray from Fabbiano, Zezas & Murray (2001), 
spectral fits from LM05. The Chandra observation was an- 
alyzed only in a central area of ~ 5.9 arcmin'^. This galaxy 
is not in the LB05 catalogue. Distance from NED (mean 
redshift-independent distance). Angular sizes and axis ra- 
tios from the RC3 catalogue; A''h from DL90. 

NGC 4144: average SFR from Ha (SFRhq = 
0.05 Mq yr^^, Moustakas & Kennicutt 2006b, using the cal- 
ibration by K98) and radio (SFRi.4GHz = O.O5M0yr"\ 
Condon et al. 2002, with the cahbration by BeU 2003);^ 
from spectra of HII regions (PT05); X-ray from LB05. 

NGC 4214: average SFR from Ha (SFRhq = 



O.18M0yr"\ Hunter 



Elmegreen 2004; SFRhc = 



O.lOM0yr ^, Schmitt et al. 2006, using the calibration 
by K98), far-infrared (SFRfir = O.15M0yr-\ Bell 2003; 
SFRfir = 0.11 M0 yr"\ Schmitt et al. 2006, using the cali- 
bration by K98) and radio (SFRi.4ghz = 0.13 M0 yr"\ Bell 
2003); Z from spectra of HII regions (PT05); X-ray from 
LB05. 

NGC 4236: SFR from radio (Condon et al. 2002, using 
the calibration by Bell 2003); Z from spectra of HII regions 
(Vigroux, Stasihska & Comte 1987); X-ray from LB05. 

NGC 4248: SFR from Moustakas & Kennicutt 
(2006b); Z from spectra of HII regions (PT05); X-ray from 
LB05. 

NGC 4254 (M 99): average SFR from Ha 
(SFRhc« = 3.9M0yr"\ Kennicutt, Tamblyn & Con- 
gdon 1994; SFRhq = 5.4M0yr-\ Buat et al. 2002; 
SFRhc« = 3.2M0yr"\ Moustakas & Kennicutt 2006b) and 
far-infrared (SFRfir = 3.48M0yr"\ Buat et al. 2002); 
metallicity gradient from spectra of HII regions (Moustakas 
& Kennicutt 2006a); X-ray from LB05. 

NGC 4258 (M 106): SFR from radio (Condon et al. 
2002, using the calibration by Bell 2003); metallicity gradi- 
ent from spectra of HII regions (PVC04); X-ray from LB05. 

NGC 4303 (M 61): average SFR from ultra-violet 
(SFRuv = 8M0yr"\ Mufioz-Mateos et al. 2007) and from 
Ha (SFRh„ = 3.6 M0 yr"\ Moustakas & Kennicutt 2006b); 
metallicity gradient from spectra of HII regions (Moustakas 
& Kennicutt 2006a); X-ray from LB05 and from Tschoke, 
Hensler & Junkes (2000), whose sources B and F are not 
present in LB05. 

NGC 4321 (M 100): SFR from Grimm et al. (2003); 
metallicity gradient from spectra of HII regions (Moustakas 



& Kennicutt 2006a); X-ray from Kaaret (2001) and from 
LB05. 

NGC 4395: SFR from ultra-violet (Mufioz-Mateos et 
al. 2007); metallicity gradient from spectra of HII regions 
(PVC04); X-ray from LM05. 

NGC 4449: average SFR from far-infrared (SFRfir = 
O.22M0yr"\ Bell 2003; SFRfir = O.21M0yr"\ RanaUi 
et al. 2003, using the calibration by K98) and radio 
(SFRi.4GHz = O.42M0yr"\ Bell 2003); Z from spectra of 
HII regions (Skillman, Kennicutt & Hodge 1989); X-ray from 
LB05. 

NGC 4485/4490 (Arp 269): average SFR from Ha 
(SFRh„ = 5.49M0yr"\ Kennicutt et al. 1987, assuming 
the calibration by K98 and an extinction of ~ 1 mag, 
Clemens, Alexander & Green 1999); metallicity gradient 
from spectra of HII regions (Pilyugin &: Thuan 2007); X- 
ray from LM05. 

NGC 4501 (M 88): SFR from radio measurements 
by Condon et al. (2002); other measurements indicate both 
lower values (1.4 M0 yr~^, from Ha photometry reported 
in Gavazzi et al. 2002), and higher values (6.0 M0 yr^^, 
from FIR measurements reported in Spinoglio et al. 2002); 
metallicity gradient from spectra of HII regions (Pilyugin et 
al. 2002); X-ray from LM05. 

NGC 4559: SFR from radio (Condon et al. 2002), us- 
ing the calibration by Bell (2003); metallicity gradient from 
spectra of HII regions (PVC04); X-ray from Cropper et al. 
(2004). 

NGC 4631 (Arp 281): SFR from far-infrared (RanaUi 
et al. 2003, using the correlation by K98); Z from spectra of 
HII regions (Pilyugin, Izotova & Sholudchenko 2008); X-ray 
from LM05. 

NGC 4651 (Arp 189): average SFR from Ha 
(SFRh„ = 1.2 M0 yr"\ Hirashita et al. 2003, using the cor- 
relation by K98), far-infrared (SFRfir = 1.6M0yr"^ BeU 
2003, using the calibration by K98) and radio (SFRi.4ghz = 
1.3M0yr^^, Bell 2003); metallicity gradient from spectra 
of HII regions (Moustakas & Kennicutt 2006a); X-ray from 
LB05. 

NGC 4656: average SFR from Ha (SFRh„ = 
O.54M0yr"\ Moustakas & Kennicutt 2006b) and radio 
(SFRi.4GHz = 1.35M0yr"\ Condon et al. 2002, using the 
calibration by Bell 2003); Z from spectra of HII regions (Pi- 
lyugin, Izotova & Sholudchenko 2008); X-ray from LB05. 

The Mice (NGC 4676, Arp 242): SFR from far- 
infrared (Brassington, Ponman & Read 2007, using the cor- 
relation by K98); X-ray estimate of Z (Read 2003); X-ray 
from Read (2003); but these observations are complete only 
for objects with Lx ^ 3 x lO^^ergs"^; so, it is quite likely 
that the actual number of ULX is larger than what we re- 
port in Table 1 (6). This galaxy is not in the LB05 catalogue. 
Distance from NED (assuming Ho — 73kms^^Mpc^^ and 
V'(LocalGroup) — 6602 kms^^). Angular sizes and axis ra- 
tios from the RC3 catalogue; A*'h from DL90. 

NGC 4736 (M 94): SFR from Grimm et al. (2003); 
metallicity gradient from Moustakas & Kennicutt (2006a); 
X-ray from LB05. 

NGC 4861 (Arp 266): SFR from Ha (Schmitt et 
al. 2006), adopting the relation in K98; Z from spectra of 
HII regions (Izotov, Thuan & Lipovetsky 1997); X-ray from 
LM05. 
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PGC 45561 (UGC 08231): SFR from Hq (Kewley 
et al. 2002), adopting the relation in K98; Z from spectra 
of HII regions (Kewley, Jansen & Geller 2005); X-ray from 
LB05. 

NGC 5033: SFR from Ha (Kennicutt & Kent 1983), 
using the correlation by K98; metallicity gradient from spec- 
tra of HII regions (PVC04); X-ray from LM05. 

NGC 5055 (M 63): SFR from radio (Bell 2003); 
metallicity gradient from spectra of HII regions (PVC04); 
X-ray from LM05. 

NGC 5194/5195 (M 51, Arp 85): SFR from ultra- 
violet (Munoz-Matcos et al. 2007) ; metallicity gradient from 
spectra of HII regions with Te method (PT05) ; X-ray from 
LM05. 

NGC 5236 (M 83): SFR from Grimm et al. (2003); 
metallicity gradient from spectra of HII regions (Pilyugin et 
al. 2006); X-ray from LM05. 

NGC 5238 (Mkn 1479): SFR from Ha (Moustakas & 
Kennicutt 2006b); Z from spectra of HII regions (Moustakas 
& Kennicutt 2006a); X-ray from LB05. 

NGC 5408: SFR from far- infrared (Stevens, Forbes & 
Norris 2002), adopting the correlation between SFR and FIR 
in K98; Z from spectra of HII regions (Masegosa, Moles & 
Campos- Aguilar 1994); X-ray from LM05. This galaxy is not 
in the LB05 catalogue. Distance from NED (mean redshift- 
independent distance). Angular size and axis ratio from the 
RC3 catalogue; Nn from DL90. 

NGC 5457 (M 101, Arp 26): average SFR from 
ultra-violet (SFRuv = 5.05 Meyr^^ Bell & Kenni- 
cutt 2001; Mufioz-Mateos et al. 2007), Ha (SFRh„ = 
2.58 Mq yr"\ Kennicutt et al. 1994; Bell & Kennicutt 2001) 
and radio (SFRi.4ghz = 2.02MQyr"\ Condon et al. 2002; 
Bell 2003; Ranalli et al. 2003) ; metallicity gradient from 
spectra of HII regions (Pilyugin & Thuan 2007) ; X-ray from 
LM05. 

Circinus (PGC 50779): SFR from Grimm et al. 
(2003); estimate of Z from optical lines (Oliva, Marconi & 
Moorwood 1999), in agreement with an X-ray estimate of 
Z (Smith & Wilson 2001). The measurement of Z might be 
affected by the central AGN; X-ray from LM05. This galaxy 
is not in the LB05 catalogue. Distance from NED (mean 
redshift-independent distance). Angular size and axis ratio 
from the RC3 catalogue; A^h from DL90. 

NGC 6946 (Arp 29): SFR from radio (Bell 2003); 
metallicity gradient from spectra of HII regions (Pilyugin et 
al. 2002); X-ray from LM05. 

PGC 68618 (IC 5201): SFR from Ha (Dopita & 
Ryder 1994; Ryder & Dopita 1994), adopting the relation 
in K98; metallicity gradient from spectra of HII regions 
(PVC04); X-ray from LB05. 

NGC 7714/7715 (Arp 284, Mkn 538): SFR from 
far-infrared (Schmitt et al. 2006), adopting the relation in 
K98; Z from spectra of HII regions (Gonzalez-Delgado et al. 
1995); X-ray from Smith, Struck & Nowak (2005). 

NGC 7742: average SFR from ultra-violet (SFRuv = 
1.00MQyr"\ Iglesias-Paramo et al. 2006), Ha (SFRh„ = 
1.40MQyr"\ Trinchieri, Fabbiano & Bandiera 1989) and 
FIR (SFRpiR = 1.42 Mq yr~\ Trinchieri et al. 1989); Z from 
spectra of HII regions (PT05); X-ray from LB05. 

Milky Way (MW): SFR from Grimm et al. (2003); 
metallicity gradient {Z{r) = 8.762 - 0.356r/i?25) from OB 
stars (Daflon & Cunha 2004); X-ray from Grimm, Gilfanov 



& Sunyaev (2002). This galaxy is not in the LB05 catalogue. 
Suggested lAU distance to MW centre. Angular size and A'h 
were not used. 

IC 10: SFR reported in table 2 of Legrand et al. (2001); 
Z from spectra of HII regions (Garnett 1990), Z from spec- 
tra of HII regions; X-ray from Wang, Whitaker & Williams 
(2005). This galaxy is not in the LB05 catalogue. Distance 
from NED (mean redshift-independent distance). Angular 
size and axis ratio from the RC3 catalogue; Ah from DL90. 

Large Magellanic Cloud (LMC): SFR from Grimm 
et al. (2003); metallicity gradient from spectra of HII regions 
(PVC04); X-ray from Long, Helfand & Grabelsky (1981). 
This galaxy is not in the LB05 catalogue. Distance from 
NED (mean redshift-independent distance). Angular size 
and Ah were not used. 

Small Magellanic Cloud (SMC): SFR from Grimm 
et al. (2003); metallicity gradient from spectra of HII regions 
(PVC04); X-ray from Wang & Wu (1992). This galaxy is not 
in the LB05 catalogue. Distance from NED (mean redshift- 
independent distance). Angular size and Ah were not used. 



APPENDIX C: ANALYSIS OF ebh 

ebh is defined as the ratio between Nulx and Nbh (see equa- 
tion 4). Table CI and Figs. CI and C2 show that the be- 
haviour of ebh as a function of the SFR and of the metal- 
licity is consistent with a constant. This result supports the 
hypothesis that Nbh is proportional to Nulx, although the 
error bars on ebh are very large. 

We point out that the trend that can be seen in the left- 
hand panels of Figs. CI and C2 (ebh versus SFR) is much 
less significant than it appears. In fact, in galaxies with a 
low SFR the presence of a single ULX can easily increase the 
value of ebh, populating the top- left part of the diagrams. 
We checked the importance of this effect by looking at the 
23 galaxies with SFR < O.3M0yr~^: the average of their 
ebh, weighted on Nbh, is ~ (9 ± 5) x 10"'*, in the P98 case. 
Such value is much more in line with the results for galaxies 
with higher SFR; it is still a bit higher than the average, but 
this might be due to some form of publication/observation 
bias. 
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Table CI. Parameters of the power-law fits and x for ebh- 



y— axis Stellar evolution model Index 



Normalization x^/dof ^ r (pi) 



SFR [Mq yr-i] 


<;bh 


P98 


SFR [Mq yr-i] 


<;bh 


P98 


Z/Zq 


<;bh 


P98 


Z/Zq 


<;bh 


P98 


SFR [Mq yr-i] 


CBH 


BIO 


SFR [Mq yr-i] 


<;bh 


BIO 


Z/Zq 


EBH 


BIO 


Z/Zq 


ffiH 


BIO 



-0.27 ±0.19 
0.00 

"•'-'^-0.27 

0.00 

-0.26 ±0.19 

0.00 

-0.16 ±0.29 

0.00 



-3 02+°-"^ 
'^•'J^-0.12 

-3.18 ±0.07 

_o -17+0.27 
•^•^'-0.23 

-3.18 ±0.07 

_o 97+0.10 
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Notes. The sample adopted for the fit in this Table is represented by all the galaxies in Table 1. ^ x^ /'^'^^ i^ the x^ divided by the 
degrees of freedom (dof). '^ r is the Pearson correlation coefficient, and pr is the probability of finding a value larger than |r|, if the two 

variables were uncorrclated and normally distributed. 
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Figure CI. ebh a-s a function of the SFR (left-hand panel) and of Z (right-hand panel) for P98. The filled circles are the galaxies listed 
in Table 1. The error bars on both the x— and the y— axis are 1 — a errors (see Section 2.3 for details). The solid lines are the power-law 
fits for the sample with Nulx > 0- The dashed lines (red on the web) are fits with a constant value for the same sample. 
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Figure C2. The same as Fig. CI for BIO. 



